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Abstract. Supersymmetric theories, including the minimal supersymmetric standard model, 
usually contain many scalar fields whose potentials are absent in the exact supersymmetric 
limit and within the renormalizable level. Since their potentials are vulnerable to the finite 
energy density of the Universe through supergravity effects, these flat directions have non¬ 
trivial dynamics in the early Universe. Recently, we have pointed out that a flat direction 
may have a positive Hubble induced mass term during inflation whereas a negative one after 
inflation. In this case, the flat direction stays at the origin of the potential during inflation 
and then obtain a large vacuum expectation value after inflation. After that, when the Hub¬ 
ble parameter decreases down to the mass of the flat direction, it starts to oscillate around the 
origin of the potential. In this paper, we investigate the dynamics of the flat direction with 
and without higher dimensional superpotentials and show that topological defects, such as 
cosmic strings and domain walls, form at the end of inflation and disappear at the beginning 
of oscillation of the flat direction. We numerically calculate their gravitational signals and 
find that the observation of gravitational signals would give us information of supersymmetric 
scale, the reheating temperature of the Universe, and higher dimensional operators. 
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1 Introduction 


The discovery of the Higgs boson [1, 2] makes it clear that the Standard Model has a hier¬ 
archy problem between the electroweak scale and the Planck scale. One of the most elegant 
solutions is a symmetry relating fermions and bosons, called a supersymmetry (SUSY), which 
results in cancellations of quadratic divergences in quantum collections of scalar field masses 
including Higgs mass. In addition, the gauge couplings are miraculously unified at the grand 
unification scale in the minimal SUSY Standard Model (MSSM). This may indicate that the 
Standard Model gauge interactions are unified as a single gauge interaction at that scale. In 
addition to these phenomenological successes in particle physics, SUSY theories have plentiful 
implications in the early Universe. 

SUSY theories, including the MSSM, usually contain many scalar fields whose potentials 
are absent in the exact SUSY limit and within the renormalizable level (see e.g., Ref. [3]). 
Since their potentials are vulnerable to the finite energy density of the Universe through 
supergravity effects, these flat directions have nontrivial dynamics in the early Universe. One 
famous example is the Affleck-Dine mechanism, in which the baryon asymmetry is generated 
through a rotational motion of a flat direction in field-value complex plane [4, 5]. In that 
mechanism, so-called a negative Hubble induced mass term plays an important role for the 
flat direction to have an initial vacuum expectation value (VEV), which is closely related 
to the amount of the resulting baryon asymmetry. Flat directions can also be considered 
as inflaton because their potentials are flat. However, inflation models in supergravity are 
usually affected by Hubble induced mass terms unless we introduce a shift symmetry [6] or 
consider a class of D-terrn inflation [7, 8] . 

In the literature, it is usually assumed that the coefficient of Hubble induced mass 
terms is constant during and after inflation [5]. However, in Ref. [9], we have pointed out a 
possibility that the coefficient of Hubble induced mass terms changes at the end of inflation. 
This is partly because in most of the SUSY inflationary models the energy density of the 
Universe is dominated by different sources between the era during inflation and the one after 
inflation. In particular, the coefficient of the Hubble induced mass term can be positive 
during inflation and negative after inflation. With such a Hubble induced mass term, the flat 
direction stays at the origin of the potential during inflation and then obtains a large VEV 
after inflation. 1 

Since the flat direction obtains a large VEV after inflation due to the negative Hubble 
induced mass term, we need to consider higher dimensional operators to stabilize its VEV. 
In some cases, the higher dimensional operator has U{ 1) or Z n symmetry. The VEV of the 
flat direction breaks the symmetry spontaneously and results in formations of topological 
defects, such as cosmic strings and domain walls, after the end of inflation. Interestingly, 
a typical width of these topological defects is always of the order of (but smaller than) the 
Hubble length because the curvature of the potential for the flat direction is of the order 
of the Hubble parameter. It follows that the energy density of the topological defects never 
dominate that of the Universe, as shown in this paper. When the Hubble parameter decreases 
down to the soft mass of the flat direction, which we expect 0(1) TeV, the potential for the 
flat direction is dominated by its soft mass term and it starts to oscillate around the origin of 

1 Note that the flat direction considered in this paper does not belong to the inflation sector. It is just one 
of the flat direction in the MSSM sector like the one assumed in the context of the Affleck-Dine baryogenesis. 
In addition, we are not restricted to generating the baryon asymmetry but consider the dynamics of general 
flat directions in SUSY theories, including the MSSM. 
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the potential. Since the symmetry is restored at the origin of the potential, the topological 
defects disappear at that time. Therefore, the topological defects exist only between the end 
of inflation and the beginning of the oscillation, so that we can not directly observe them 
at the present epoch. However, we have a chance to observe redshifted gravitational waves 
(GWs) emitted by the topological defects. 

Stochastic GW signals are generated from the dynamics of cosmic strings [10-13] and 
domain walls [14-16]. In the present case, emission of GWs begins at the end of inflation 
and ceases when the Hubble parameter decreases down to the soft mass of the flat direction. 
Since the emitted GW signals have a peak at the wavenumber corresponding to the Hubble 
length, resulting GW signals have a peak at the wavenumber corresponding to the soft mass 
of the flat direction. This implies that in principle we can obtain the information of the 
SUSY scale from the observation of GW signals [9]. In addition, since the energy density of 
the topological defects is roughly proportional to the square of the VEV of the flat direction, 
which is related to its higher dimensional operators, its GW signal also has information of 
higher dimensional operators. Finally, since the GW spectrum is sensitive to what dominates 
the energy density of the Universe, we could obtain information of the reheating temperature 
from the observation of GW signals [17, 18]. 

In this paper, we investigate the dynamics of the flat direction which obtains a positive 
Hubble induced mass term during inflation and obtains a negative one after inflation. We 
consider three cases for its higher dimensional operators to stabilize its VEV; the one coming 
from a nonrenormalizable superpotential, the one coming from a Kahler potential with 17(1) 
symmetry, and the one coming from Kahler potentials with Z n symmetry. We numerically 
follow its dynamics and confirm the above qualitative understanding, i.e., topological de¬ 
fects, such as cosmic strings and domain walls, form after the end of inflation and disappear 
when the Hubble parameter decreases down to the soft mass of the flat direction. From 
the simulations, we also evaluate its GW signals and discuss its detectability by future GW 
detectors. 

In the next section, we briefly review flat directions in SUSY theories. In Sec. 3, we 
explain the potentials for flat directions originating from SUSY breaking and higher dimen¬ 
sional operators. We consider three important cases for higher dimensional operators and 
show our results of numerical simulations for each case in Secs. 4, 5, and 6. Finally, Sec. 7 is 
devoted to the conclusion. 


2 Flat directions in the MSSM 


A large number of scalar fields are introduced in SUSY theories because SUSY relates 
fermions and bosons. The potentials for scalar fields are severely restricted by SUSY, which 
especially results in existence of many flat directions. Flat directions are scalar fields whose 
potentials are absent within the renormalizable level as long as SUSY is unbroken. 

Here, we illustrates how the potentials for scalar fields are absent in SUSY theories by 
taking a flat direction called the u c d c d c flat direction as an example. Let us focus on a scalar 
field constructed by right-handed squarks through the following orthogonal matrix: 





( 2 . 1 ) 


where the lower and upper indices represent flavor and color, respectively. The dots represent 
other directions, which we are not interested in. Since the inverse matrix is given by the 
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transposed matrix, we obtain 





In the MSSM, the superpotential is given by 2 


Wmssm = y u H u Qu c - y c iH d Qd c - y e H d Le c + yH u H d , 


( 2 . 2 ) 


(2.3) 


within the renormalizable level. Here we have omitted flavor indices for notational simplicity. 
It is easy to see that the F-term potential for the field (p is absent such as 


V F = 


SHmssm 


dtp 


dl^MSSM 


d(u c )i 


R 


+ 


cWmSSM 


d{d c ft 


G 


+ 


^IHmssm 


d(d c 


= 0. 


The F-term potential is also absent like 


\D%\ 2 = gl | {u c )i*T a (u c )i + ( d c )f*T a (d c )f + (d c )f*T a (d c )f | 2 


= ^|^>| 4 |Tr (T a )| 2 = 0, 


l^il 2 = 


-| |(^ C )f I' + | |(^ C )f | + | |(^ C )f | 


= 0, 


(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 


where |-D!j| 2 and \Di\ 2 are H-term potentials for SU(3 ) and U(l)y, respectively. Therefore, 
the field cp has a flat potential and is called a flat direction. The above example consists of 
right handed squarks of u c ,d c ,d c , and is called the u c d c d c flat direction. It is known that 
every flat direction is characterized by gauge-invariant monomial in this manner. Table 1 is 
the list of flat directions in the MSSM [3], Note that there are many flat directions even in 
such a simple model. 

Although flat directions have no potential within the exact SUSY limit and the renor¬ 
malizable level, they obtain nonzero potentials through SUSY breaking and nonrenormaliz- 
able operators (i.e., underlying higher energy theory). These potentials induce non-trivial 
dynamics of flat directions. The next section is devoted to discussing this point. 


3 Potentials for flat directions 

In this section, we discuss the induced potentials for flat directions through SUSY breaking 
and nonrenormalizable operators. Flat directions have soft (SUSY breaking) masses of the 
order of sparticle masses, which are subject to collider experiments and should be larger than 
0(1O 2 ' 3 ) GeV [19, 20]. In addition, the finite energy density of the Universe contributes to 
potentials for flat directions. For instance, scalar fields obtain so-called Hubble induced terms 
through supergravity effects during inflation because inflation is driven by a finite energy 
density [5]. This is also the case during the inflaton oscillation dominated era as we explain 
in Sec. 3.2. Finally, when a flat direction has a large VEV, nonrenormalizable operators 
become important. 3 In this paper, we consider three important cases for nonrenormalizable 
operators, which we specify in Sec. 3.3. 

2 We assume a symmetry, such as R-parity, to forbid baryon and lepton number violating renormalizable 
terms. 

3 Although thermal effects sometimes play an important role to investigate the dynamics of flat directions [5, 
21-23], it is irrelevant in this paper (see footnote 8). 
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Table 1. Flat directions in the MSSM and possible values of the power of higher dimensional 
superpotentials [3]. 


flat directions n 


LH U 

4, 

6, 8, ... 

H u H d 

4, 

6, 8, ... 

u c d c d c 

6, 

9, 12, ... 

LLe c 

6, 

9, 12, ... 

Qd c L 

6, 

9, 12, ... 

QQQL 

4, 

8, 12, ... 

Qu c Qd c 

4, 

8, 12, ... 

Qu c Le c 

4, 

8, 12, ... 

u c u c d c e c 

4, 

8, 12, ... 

d c d c d c LL 

5, 

10, 15, ... 

u c u c u c e c e c 

5, 

10, 15, ... 

Qu c Qu c e c 

5, 

10, 15, ... 

QQQQu c 

5, 

10, 15, ... 

{QQQ) 4 LLLe c 

7, 

14, 21, ... 

u c u c d c Qd c Qd c 

7, 

14, 21, ... 


3.1 Soft terms 

Since we have not discovered SUSY particles yet, SUSY particles have to be much heavier 
than the SM particles. This can be achieved by introducing a SUSY breaking hidden sector, 
which induces soft mass terms to SUSY particles through some mediation mechanism. The 
potential for a flat direction can be written like 

V{ct>)=ml H 2 , (3.1) 

around the weak scale vacuum. Note that not only soft masses but //-term can also contribute 
to rrifp if the flat direction consists of Higgs field. 1 Here we neglect A terms because they are 
subdominant and irrelevant in the following discussion. 4 5 In this paper, we assume that weak 
scale SUSY is broken at latest earlier than the time at which the Hubble parameter H(t ) 
decreases down to so that the soft mass term of Eq. (3.1) dominates the potential for the 
flat direction after H(t) ~ m^. 

3.2 Hubble induced terms 

The potential for a flat direction cj) is modified by supergravity effects during and after 
inflation. The energy density of the Universe is so large that we have to consider its effects 

4 In this paper, hereafter, we neglect the flavor indices. This is valid when the Hubble induced terms and/or 
nonrenormalizable terms breaks flavor symmetry completely. The case with a flavor symmetry is considered 
in Appendix C. 

5 There might be Hubble induced A terms during inflation, which plays important role in the context of 
the Affleck-Dine baryogenesis. In the present case, the flat direction stays at the origin during inflation, so 
that such higher dimensional terms are irrelevant. In addition, Hubble induced A terms are irrelevant after 
the end of inflation because the time-averaged F term (F) vanishes (see e.g., Ref. [24]). Therefore, in the 
present scenario, the flat direction can not generate the baryon asymmetry. 
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on the potential for the flat direction. In the supergravity, scalar potentials are written in 
terms of superpotential, W, and Kahler potential, K , which are functions of scalar fields i/j 1 
specified below. The potential for scalar fields is given by 


V = e K ! M pi 


DiWK ij DjW* 


3\W\ 2 /M^ , 


(3.2) 


where K 1 ^ = ( K jj) -1 = (c^K/ch/^ch/;*^) -1 and DjW = dW/d^+W/Mp^dK/d^ 1 . Hereafter, 
we assume E-term inflation models (see Appendix A for example), in which we can neglect 
the contribution from D-term potentials. The Kahler potential also determines kinetic terms 
such as 


£ kinetic = K fj 8^(8^)*. (3.3) 

One may assume a minimal Kahler potential of K = | -0* | 2 , which results in the canonical 
kinetic term. However, since we do not have any knowledge about the physics above a cutoff 
scale, such as the Planck scale, there are expected to be cutoff suppressed terms in the Kahler 
potential (see Eq. (3.4)). 

The relevant fields are a flat direction 0 and two fields in the inflaton sector denoted 
as / and X, so that 0* are identified as 0, /, and X in Eqs. (3.2) and (3.3). The superfield 
I contains a field x which starts to oscillate after inflation and whose energy dominates 
the energy density of the Universe, that is, % is identified as inflaton in chaotic inflation 
models [6] and as water-fall fields in hybrid inflation models [25, 26]. The minimal Kahler 
potential for the field I is given by K rn \ n (f) = (/ + 1 *) 2 /2 for the case of the chaotic inflation 
model proposed in Ref. [6] and K m i n (I) = \I\ 2 for the case of hybrid inflation models (see 
Appendix A). The superfield X represents the field whose E-term potential energy drives 
inflation. In other words, the E term of X satisfies the relation of \dW- m {/dX] 2 = 2>H 2 nf Mpp 
In Appendix A, we briefly explain chaotic and hybrid inflation models and identify the fields 
I and X in those models. 

To see how the potential for the flat direction is modified, let us consider the following 
non-minimal Kahelr potential: 

K = |0| 2 + |A| 2 + || |A| 2 H 2 + K min (J) (l + ^ H 2 ) , (3-4) 

where M* is a cutoff scale, and c a (a = X , I) are constants of order one. We expect that M* 
is less than or equal to the Planck scale Mp\ (~ 2.4 x 10 18 GeV). We set M* = Mpi in Sec. 4, 
while we assume M* to be less than Mp\ in Secs. 5 and 6. Since the F term of X satisfies 
the relation of \dW/dX\ 2 = 3H 2 , f Mp l5 the flat direction obtains the following mass, called 
the Hubble induced mass, during inflation: 

= c Hint H 2 ni H 2 , (3.5) 

°H in f = 3 ^1 - cxjff'j ■ (3.6) 

Let us note that cn in{ is of the order of MpjM 2 when cx is of order one. While the coefficient 
CH inf is assumed to be negative in the context of the Affleck-Dine baryogenesis [4, 5], we 
assume it to be positive in this paper. In this case, the flat direction stays at the origin, i.e., 
0 = 0, during inflation. 
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After inflation ends and before reheating completes, the energy density of the Universe is 
dominated by that of the oscillation of a scalar field (denoted by x) and the Hubble parameter 
decreases with time as H(t) oc a _3 / 2 (f), where a(t) is the scale factor. The field x oscillates 
with the frequency of its mass m x . The oscillation amplitude is redshifted as 


X oc Xi 



(3.7) 


where Xi is the oscillation amplitude just after the oscillation starts around t = t{. Since the 
oscillation time scale of x is much smaller than H _1 , that is, H <C rn x . we can average \' 2 
and x 2 over the oscillation time scale like 

<X 2 } ^ {m 2 x x 2 ) - (\dW/dX | 2 ) ~ \h 2 M 2 v (3.8) 

The flat direction obtains a Hubble induced mass through higher dimensional kinetic terms 
as well as the E-term potential of Eq. (3.2) during this oscillation dominated era. From 
Eq. (3.3), we can see that the relevant term comes from 

K nW , 2 d ^ 2 H 2 ’ ( 3 - 9 ) 


which, together with the one coming from the E-term potentials, results in the Hubble 
induced mass term of cuH 2 (t ) \<fi\ 2 with 


c H nsa = - 1 -c x 


M 2 


3 M 2 

2 C/ W’ 


(3.10) 


during the oscillation dominated era. Let us note that c# osc is of the order of Mpj/M 2 when 
cx and c/ are of order one. 

As we can see from Eqs. (3.6) and (3.10), the coefficient of the Hubble induced mass term 
during inflation, Cff inf , is generally different from the one during the oscillation dominated 
era, ch osc [9]. This is because the held A, whose E term drives inflation, is generally different 
from the held /, whose oscillation energy dominates the energy density of the Universe after 
inflation. In Appendix A, by taking the models of chaotic inflation and hybrid inflation in 
supergravity as examples, we see that the held X is indeed different from the held I. While 
the nonrenormalizable terms are heavily dependent on the high energy physics beyond the 
cutoff scale M*, we assume that among many hat directions in SUSY theories (at least) one 
hat direction has cn ini > 0 and c# osc < 0. In this case, the hat direction stays at the origin 
(cj> = 0) during inflation, and then, it starts to roll down to a large VEV due to the negative 
curvature after inflation. 

Here we comment on how natural to consider the above hat direction. There is no 
guiding principle to determine the signs of Hubble induced mass terms unless we specify the 
underlying high energy model beyond cutoff scale. This is what the well-known Affleck-Dine 
mechanism is also based on, while it assumes a negative coefficient during inflation. As shown 
above, the coefficient during inflation and that after inflation do not have the common origin 
and hence can be different in general. We consider one out of four possible combinations of 
the signs of the two coefficients. It is natural that one or more hat directions have such a 
combination because there are a lot of hat directions in the MSSM (see Table 1). This is 
why our work can find its application in not all but wide class of supersymmetric models. 
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3.3 Higher dimensional terms 

Since the flat direction obtains a large VEV due to the negative Hubble induced mass term 
after inflation, we need to consider its higher dimensional terms to stabilize the VEV of the 
flat direction. Depending on symmetries that the flat direction possesses, three important 
cases are considered in this paper. 

In the next section, we consider the case that the flat direction has a nonrenormalizable 
superpotential like 


W = A 


nM i 


n —3 ’ 

PI 


(3.11) 


where n (> 4) is an integer dependent on each flat direction and A is an positive integer. 
Let us stress that among possible higher dimensional terms, only the lowest dimensional one 
is relevant to the dynamics of the flat direction. In the case of the u c d c d c flat direction, for 
instance, it may come from a term like (see Eq. (2.2)) 


W = A 


(i u c d c d c f 

6Mp, 


A 


</> 6 

6 M|, 


(3.12) 


In this case, the power of superpotential n is 6. Note that n = 3 superpotential is forbid¬ 
den due to R-parity. When we extend this symmetry to a discrete R-symmetry, n can be 
6, 9,12,... for the u c d c d c flat direction. For each flat direction, possible values of n are shown 
in Table 1 [3]. The E-term potential for the flat direction comes from the superpotential and 
is given by (see Eq. (3.2)) 


Vnr - 


dW 

dcj) 


\2n —2 




(3.13) 


where we neglect irrelevant higher dimensional terms. Note that the potential of Eq. (3.13) 
possesses a global 17(1) symmetry. This is related to the baryon charge (or the lepton charge 
or a combination between them), which plays an important role in the context of the Affleck- 
Dine baryogenesis. Let us emphasize that the symmetry is a global symmetry independent 
from the local symmetries of the flat direction. The symmetry can be explicitly described as 
(see Eq. (2.2)) 


( (u c ) f —> e l6 {u c ) f 

(/> —»• e i0 cj) ^ l (d c )f -> e ie (d c )f , (3.14) 

1 ( d c )f -> e i9 (d c ) f 

for the case of the u c d c d c flat direction. One can see that this transformation can not be 
described by the gauge symmetries of 517( 3) x 17(l)y, so that the 17(1) symmetry is an 
additional global symmetry. 

If A in Eq. (3.11) is sufficiently small or superpotential is absent due to some symmetry, 
the VEV of the flat direction can be as large as the cutoff scale M*. For example, if we intro¬ 
duce B — L symmetry, it forbids superpotentials for the u c d c d c flat direction like Eq. (3.12). 
If we assign R-charges as, for example, R(Q) = R(L) = R(e c ) = R(u c ) = R(d c ) = 1 and 
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R(H U ) = R(Hd) = 0, it also forbids superpotentials for the u c d c d c flat direction. 6 In such 
cases, the potential for the flat direction may be dominated by higher dimensional terms 
coming from nonrenormalizable Kahler potential [27]. The induced potential is proportional 
to H 2 like the Hubble induced mass term described above. In Sec. 5, we consider a flat 
direction with a potential given by 


\i\zm—z 

H N r(</>) = a H H 2 (t y ' 4 , 

iVipi 

where m (> 3) is a certain integer. ' The coefficient an is given by 


(3.15) 


=a ~U:J 


2m—2 


a H — &H 


(3.16) 


where we assume a’ H is positive and of order one. Note that the potential of Eq. (3.15) 
respects 17(1) symmetry. 

On the other hand, if we assign R-charges as, for example, R(Q) = R(L) = R(e c ) = 
R(u c ) = R(d c ) = 0 and R(H U ) = R(Hd) = 2, nonrenormalizable Kahler potentials may not 
respect 17(1) symmetry. In Sec. 6, we consider the same potential with Eq. (3.15) except for 
an additional 17(1) breaking term coming from higher dimensional Kahler potentials: 


12m—2 


Vnr(0) = a H H 2 (t) — f2m _ 4 - 
pi 


Mi 


m'—2 

PI 


M n 


+ c.c. 


(3.17) 


where m! = 2, 3,4,..., is an integer. We assume m! < 2m — 2 so that the potential for the 
flat direction is bounded from below. The coefficient bn is given by 


b H 


b\ 


H 



(3.18) 


where we assume b' H = 0(1). Here we have redefined the field <j) such that b' H is postitive. 
This additional term explicitly breaks the 17(1) symmetry, but there remains Z m i symmetry 
in that potential. 


3.4 Summary of this section 

Here we summarize the potential for the flat direction. It is given by 

V(<j>) = ml H 2 + c H H 2 (t) \ct>\ 2 + K N r, (3.19) 

c The R-symmetry has to be broken to obtain nonzero gaugino and higgsino masses and to set the vacuum 
energy (almost) zero. Its symmetry breaking order parameter is proportional to the SUSY breaking scale, so 
that R-symmetry breaking terms are suppressed by a factor of m|/M* . These terms may be comparable to 
the other terms, such as the soft mass term, around the time when the Hubble parameter decreases down to 
the soft mass scale. This implies that there is 0(1) uncertainty in our result when we omit those R-symmetry 
breaking terms. On the other hand, it might be possible that R-symmetry breaking terms are suppressed 
by a factor of m^,/Mp\. In this case, one can easily find that they can be neglected even when the Hubble 
parameter decreases down to the soft mass scale. After that, higher dimensional terms are irrelevant because 
the potential for the flat direction is dominated by the soft mass term. Therefore, it is reasonable to consider 
that R-symmetry may suppress the superpotential for the flat direction. 

7 In general, two or more terms can give comparable contributions to the potential for flat direction. While 
it can change our results by a factor of (7(1), in this paper, we consider a nonrenormalizable potential with a 
single term for simplicity. 
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where the higher dimensional terms Vnr is given by Eqs. (3.13), (3.15), or (3.17). Hereafter, 
Cff inf and ch osc are collectively denoted as ch■ We consider a flat direction which has a 
positive Hubble induced mass term (c# = c# inf > 0) during inflation and a negative one 
(cjj = ch osc < 0) after inflation. 

In the following sections, we investigate the dynamics of the flat direction and show that 
topological defects, such as cosmic strings and domain walls, form after the end of inflation. 
We numerically follow the dynamics of the flat direction and calculate its GW signal. We also 
discuss the detectability of GW signals and explain that the information of the mass of the 
flat direction, the reheating temperature of the Universe, and higher dimensional potentials 
are imprinted on GW signals. In Secs. 4, 5, and 6, we consider the case that the flat direction 
has a higher dimensional potential of Eqs. (3.13), (3.15), and (3.17), respectively. 


4 Case A: W ^ 0 


In this section, we consider the case that the flat direction has a nonrenormalizable superpo¬ 
tential of Eq. (3.11). The E-term potential for the flat direction is given by 


Cnr 


dW 

d<j> 


\2n —2 


M, 


2n—6 ’ 
PI 


(4.1) 


where n > 4. 

In the next subsection, we briefly explain the dynamics of the flat direction in the early 
Universe. We numerically follow its dynamics and calculate its GW signal, and the results 
are shown in Sec. 4.2. Then we give physical interpretation of the results in Sec. 4.3. Finally, 
we discuss the detectability of GW signals in Sec. 4.4. 


4.1 Dynamics of the flat direction 


The equation of motion for the flat direction is given by 

4> + 3H<j ) -\v 2 0+^ = O, 
a z op* 


(4.2) 


with dots denoting derivatives with respect to t. The potential for the flat direction V((j>) is 
written as 

I /12m—2 

V&) = ml H 2 + c H H\t ) H 2 + (4.3) 

Mpi 

Note that we collectively denote c# inf and ch osc as ch ■ As explained in the previous section, 
we consider a flat direction which has a positive Hubble induced mass term (c# inf > 0) during 
inflation and a negative one {ch osc < 0) after inflation. During inflation, the flat direction 
stays at the origin, = 0, due to the positive Hubble induced mass term. Since the coefficient 
ch becomes negative after the end of inflation, the flat direction rolls down to the stable point 
and becomes to obtain a large YEV. 8 The VEV is given by 


\ch osc \ 

A 2 (n — 1) M 2 ! ) 


1 /( 271 - 4 ) 


Ml 


Pll 


(4.4) 


8 One might wonder if the flat direction is affected by thermal effects and is forced to stay at the origin 
by the effective thermal mass for the flat direction [5, 21-23]. However, we can neglect thermal effects on the 
flat direction by the following discussion. During the inflaton oscillation dominated era, the temperature of 
the plasma is given by T ~ (H(t)rjMpi) 1 ' 4 , where Tj is the decay rate of the inflaton. The flat direction 
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as long as \c Hosc \H 2 (t) > m|. 

The potential of Eq. (4.1) has a global U( 1) symmetry, which corresponds to the baryon 
charge (or the lepton charge or a combination of them) in the context of the Affleck-Dine 
baryogenesis. After the end of inflation, the non-zero VEV of the flat direction breaks the 
global U( 1) symmetry. 9 Since the flat direction stays at the origin during inflation, cosmic 
string network forms after this phase transition. The energy of cosmic strings per unit length 
p is roughly given by 

M~(I0I> 2 , (4-5) 

depending on dimensionless parameters only logarithmically [28]. By using numerical sim¬ 
ulations (see the next subsection), we confirm that cosmic string network reaches a scaling 
regime well before the string disappears. In this regime, the number of cosmic strings in the 
Hubble volume is 0(1) and the energy density of cosmic strings p cs is given by 

H~ l x H 3 = pH 2 . (4.6) 

A typical width of cosmic strings is roughly given by the curvature of the potential, which is 

of the order of the Hubble length in the present case. 10 This means that the width of cosmic 
strings increases with time. This is one of the outstanding characteristics of cosmic strings 
in this scenario compared with other cosmic strings considered in the literature. 

Since the Hubble induced mass squared decreases with time as oc a~ 3 (t) after the 
inflation ends, the soft mass term rn 2 . \cj)\ 2 eventually dominates the potential for the flat 
direction and the flat direction starts to oscillate around cj> = 0. This implies that cosmic 
strings disappear at the time of |cfl- osc | H 2 (t) — m which can be rewritten as 

t - tdecay = 2 ^ CHosc ^ . (4.7) 

3 

With 7B0 = 0(1) TeV, cosmic strings disappear much earlier than the bing bang nucle¬ 
osynthesis (BBN) epoch, so that observational constraints on long-lived topological defects 
are not applicable to our scenario. The flat direction decays into radiation through gauge 
interactions soon after cosmic strings disappear. 

We define the reheating temperature Trh as the temperature when the energy density 
of the Universe is dominated by that of radiation. The Hubble parameter at that time is 
given by 

< 4 - 8 > 

obtains a thermal mass of the order of the temperature. The ratio of the thermal mass to the Hubble induced 
mass is thus roughly given by T/H(t) ~ (YiM^i/H(t) 3 ) 12 ^■ This is smaller than unity just after inflation if 
we consider a Planck suppressed decay of inflaton. Thus, the flat direction obtains a large VEV soon after 
the end of inflation due to the negative Hubble induced mass term. Once the flat direction obtains a large 
VEV, thermal effects become suppressed and can be neglected because particles in thermal bass obtain large 
effective masses [22, 23]. Therefore, we can neglect thermal effects on the flat direction. 

9 It also breaks SM gauge symmetries because flat directions have nonzero SM gauge charges. In the 
unitary gauge, only gauge fields may obtain a nontrivial configuration. Since the broken gauge field obtains 
the effective mass of g {(/>), the gauge field configuration begins to oscillate and becomes trivial = 0) soon 
after the phase transition. Therefore, we focus only on the global symmetry of the flat direction. See also the 
discussion in Appendix C.l. 

10 Since a typical width of cosmic strings is of the order of the Hubble radius, the Nambu-Goto approxima¬ 
tion, which was used in Ref. [29] for example, is inappropriate to describe these cosmic strings. 
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where g * is the effective relativistic degrees of freedom for the energy density. Note that in low- 
scale SUSY models we should require Trh 1$ 10 9 GeV to avoid the gravitino problem [30, 31]. 
Even for high-scale SUSY models such as pure gravity mediation, Trh ;$ 10 10 GeV is required 
to avoid an overproduction of LSP (= 0(100) GeV) [32-34], In such well-motivated cases, the 
cosmic strings disappear (iLdecay = (2/3)tj e 1 cay = 0(1O 2 ' 3 ) GeV) before reheating completes 
(77 rh 1$ 0(1) GeV). Hereafter we consider such a case. 

After they form at the end of inflation and before they disappear at t ~ fdecay) cosmic 
strings emit GWs. We numerically follow the evolution of the flat direction and calculated 
its GW signal. The results are shown in the next subsection. 


4.2 Results of numerical simulations 

In this subsection, we show results of our numerical simulations for the dynamics of the 
flat direction and its GW signal. Since we consider the evolution of cosmic strings during 
the oscillation dominated era (H <C m x ), the Hubble parameter decreases with time as 
H oc a -3 / 2 . 

By redefining the variables, we rewrite the equation of motion such as 




// 


dx 2 






4 

+ (n- 



|^-V = 0,(4.9) 


where 


ip = (A -1 Mpf 3 Rj)~ 1/(n-2) a(t)</> 

= (\CH 0 J /{n - l)) 1 /(2n-4) (l&l) , (4.10) 

dr = H{dt/a(t), (4-11) 

m<j> = m^/Hi (< 1), (4.12) 

x = HiX. (4.13) 


Let us emphasize that the parameter A can be absorbed by the redefinition of (f>, so that our 
numerical simulation can be applied to any cases with an arbitrary value of A. The variable 
t is a dimensionless conformal time and is related with comoving horizon size r through 
r = TdiHi . Primes denote derivatives with respect to r. We rewrite important parameters 
by the conformal time and obtain 


n = 2 , 


H(t) = ( | ) Hi , 


7~decay 


|cf/oJ H, 


2 \V 6 


Ti , 


771 : 


^RH = 


/ 90Mp,iL| 

v 


2 \ 1/6 


RH 


Ti, 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


where we implicitly assume a* = a(r)) = 1. The subscript i represents the values at the time 
of 7 at which we set initial conditions for numerical simulations. 

We have performed three-dimensional lattice simulations to solve the classical equation 
of motion Eq. (4.9) numerically. We use a numerical method similar to the one used in 
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Ref. [16]. They use the 4-th order symplectic integration method to evolve the equation 
of motion for the flat direction in the expanding background, while we have used the 6-th 
order symplectic integration method. We ignore the backreaction of metric perturbations 
on the field evolution, which is negligible for p cs H 2 (t)Mp V that is, for (0) Mpj (see 
Eqs. (4.5) and (4.6)). In the set of simulations, we take a unit of Hi = H(t,) = 1 and set 
m^/Hi = 5 x 10 -4 and Cff osc = —15. The final time ry and the time step At is set to be 50 
(100) and 0.2 (0.1), for N = 128 3 (256 3 ) grid points, respectively. The comoving box size is 
L = 100 (200) with N = 128 3 (256 3 ) grid points, which implies that the comoving grid size is 
Ax = 100/128(= 200/256) ~ 0.75. Since a typical width of cosmic strings is of the order of 
the Hubble radius 7/ _1 (r) = (r/r*) 3 , it is correctly resolved in our numerical simulation. We 
use the periodicity boundary condition for spatial directions. We require that the horizon 
scale is smaller than the simulation size (i.e., L/2 > r) to avoid the boundary effect on the 
cosmic string network. The initial condition for simulations is seeded by quantum fluctuations 
around the origin of the configuration space. To be concrete, we set initial conditions in the 
following way. First, we draw two samples /- i, ff: 2 f rom a normal distribution with variance 
of 


( 



(kH„.l/("-i)) 1 /( ”~ 2 ) Jv a 

2 kV> (|0i|) 2 


(4.18) 


for each discretized wavenumber vector k = (2ir/L) n, where n is a triplet of integers from 
—N/2 to N/2 — 1. In the above expression, we set A such that (|0j|) = 10 2 in units of 
Hi. Let us stress that once the cosmic string network follows the scaling evolution, the 
initial condition (i.e., A) becomes irrelevant. We generate initial conditions by taking 

a superposition of (discrete) Fourier components given by 

= 7! (4 i + 4 2 ) * ( 4 - 19 ) 

^(k) = (/*,! - 4a) • (4.20) 

Figure 1 shows one example of cosmic string network evolution obtained by our numer¬ 
ical simulation. The blue region represents the inside of cosmic strings, which we define by 
|</>(x)| < (|0|) /5. One can find that a typical width of cosmic strings increases with time 
and is roughly proportional to the Hubble length. The black line at the bottom-left corner 
in Fig. 1 describes a unit of the comoving horizon size (given by t) at each time. One can 
see that each Hubble volume contains 0(1) cosmic strings, which means that cosmic string 
network is in a scaling regime. From our numerical simulation, we find that the spatially 
averaged magnitude of 0 grows with time and eventually reaches of the order of the VEV 
(Eq. (4.4)) around Tf orm ~ 10 for the case of ch osc = —15. 

From the numerical simulation for cosmic string dynamics, we can obtain the transverse- 
traceless part of the anisotropic stress by the relation of 

x) = Ajj jim (x) (<9;0*<9 m 0 + c.c.), (4.21) 

where the projection operator in Fourier space is given by 

A ij,im(k) = Pu(k)P jm (k ) - ^ Pij(k)Pi m (k ), (4.22) 

Pij(k) = dij - kjkj, (4.23) 
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Figure 1. Cosmic string networks obtained by our numerical simulation at TGq-ff; = 20 (upper left 
panel), 50 (upper right panel), 80 (lower left panel), 100 (lower right panel). We assume ch osc = —15 
and n = 6. In the blue regions, the magnitude of the flat direction satisfies |</>(x)| < (|0|) /5, where 
(|</>|) is the VEV given by Eq. (4.4). The black line at the bottom-left corner represents the length of 
the comoving horizon size r at each time. 


with k = \k\k and Kronecker’s delta 6ij. Given trh, we can calculate the energy density of 
GWs from Eqs. (B.6), (B.7), and (B.14) in Appendix B. We define the energy spectrum of 
GWs p gw such as 


figw(r) 


1 dp gw (r) 
Ptot (t) d log k 


(4.24) 


where ptotCr) (= 3Mpji4 2 (r)) is the total energy density of the Universe. Figure 2 shows 
evolution of GW spectra obtained from our numerical simulations. The GW energy density 
grows as a typical value of the flat direction increases until it reaches the potential minimum 
at r = Tf orm ~ 10. Then the GW peak energy density decreases with time as f2 gw oc r - 12 /( n - 2 ) 
and the peak wavenumber fc pea k decreases with time as fc pea k oc r' 1 for the case of n > 8. 
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Figure 2. GW spectra obtained from our numerical simulations. The energy density of GWs are 
rescaled by (Hi/ AM P i) _4 A n_2 ). The red, green, and magenta curves are GW spectra at the conformal 
time of TCLiHi = 10, 20, and 50, respectively. We also show the time evolution of the GW peak location 
as the blue solid (dashed) line for the case of /H, = 0 (5 x 10 -4 ). We take n = 6 (upper left panel), 
8 (upper right panel), 10 (lower left panel), and 12 (lower right panel), fixing ch obc = —15. The black 
curves are plotted by Eq. (4.33) at raiHi = 50 and well describes large wavenumber modes. The gray 
dotted curves are estimated from the discussion of causality at raiHi = 50. 


On the other hand, the GW peak energy density decreases with time as n gw oc r~ 2 and its 
peak wavenumber is constant in time for the case of n = 6. We give physical interpretation 
of these results in the next subsection. 

Taking a non-zero value of into account, we find that the flat direction starts to 
oscillate around the origin of the potential and cosmic strings disappear eventually around 
t = Tfjecay — 40, which is given by Eq. (4.16) with rri^ / H t = 5 x 10~ 4 and c# OBC = —15. After 
that, the GW spectrum is freezed out, which means that the GW peak wavenumber becomes 
constant and the GW spectrum decrease adiabatically as oc r” 2 (oc a -1 ). Note that we have 
taken = 5 x 10 -4 above just for computational time saving, though it does not change 
the dynamics of flat direction until r ~ Td ecay . Since Hi is roughly given by m x Xi/Mp\ (see 
Eq. (3.8)), m<f) can be much smaller in general. When we take rn x = 10 13 GeV and Xi = Afp] 
motivated by chaotic inflation model and m^/ |cff osc | 1//2 = 10 3 GeV, Td eC ay — 4000 is larger 
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than Tf orm ~ 10 by three orders of magnitude. 


4.3 Physical interpretations 

In this subsection, we discuss how we can understand the results of our numerical simulations 
shown in Fig. 2. 

During the scaling regime, the typical length scale of cosmic string dynamics is the 
Hubble length. As a result, the emission peak wavenumber of GWs emitted from these 
cosmic strings is given by the Hubble scale: 

k A 

~ H(t). (4.25) 

a(t) 

We can estimate the energy density of GWs in the following way. Its quadrupole moment Q 
can be roughly given by [35] 

Q ~ H~ 2 x fiH- 1 , (4.26) 


where [iH 1 is the total energy of cosmic strings within a Hubble volume. Such cosmic 
strings emit GWs with the luminosity of 

C ~ M~ 2 Q 2 ~ M~ 2 (H 3 Q) 2 , (4.27) 


where we replace the time derivative with H in the last line. Thus, the produced energy 
density of GWs can be estimated as 



H 3 x H~ l JZ 


A logr 

peak H 2 M 2 , 

V M P] ) 


I I ir2w\ V(® 2) 

I c H osc \ H {t) \ 

A 2 (n — 1) Mpj ) 


(4.28) 

(4.29) 


where we use Eqs. (4.5) and (4.4) in the second and last equality, respectively. During each 
Hubble time, cosmic strings emit GWs with the emission peak wavenumber of Eq. (4.25) and 
with the produced energy density of Eq. (4.29). 

We have measured ADgw/Alogr from our lattice simulation as shown in Fig. 3. We 
determine the numerical prefactors for the emission peak wavenumber and produced energy 
density of GWs from the numerical results such as 


k A 

peak 

a(f) 


~ 2.5 H(t), 


AH 


gw 


Alogr 


10 


- peak 


C H 


H- 


'm 


2/(n—2) 


A 2 (n — 1) M 2 ! J 


(4.30) 

(4.31) 


We find that the numerical prefactors are independent of n and r (see Fig. 3), so that 
Eqs. (4.30) and (4.31) can be used for any n and r. 

From Fig. 3, we can find that AH gw /Alogr is proportional to a certain power of k for 
large and small wavenumber modes: 


Alogr 



A log T 

peak ^ 


gw 

A log r 

peak ^ 


for k < k A ak 

for k > k A ak 


(4.32) 
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Log [k /( a H (t))] 




Log[& /(a H (t))] 



Figure 3. Derivatives for GW spectra in terms of conformal time. The energy density of GWs are 
rescaled by the time dependent factor given by Eq. (4.29). We plot GW spectra at the conformal 
time of TdiHi = 20,30,40,50,60,80, and 100 for n = 6 (upper left panel), 8 (upper right panel), 10 
(lower left panel), and 12 (lower right panel), fixing ch obc = —15. We fit small and large wavenumber 
modes by the dot and dashed red lines, which are proportional to k 2 and fc -2 , respectively. 


where a ~ 2. We also find that these values are almost independent of n and r from Fig. 3. 
The GW spectra for large wavenumber modes in Fig. 2 is expected to be given by integrating 
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the GW peak energy density from Tf orm to r, such as 


^gw ( t ) 


d log T 


,a(r') Afi gw (r / ) 


' Uorm 

r^k 

'norm 


d log T 


+ 


d log r' 


a(r) 

A log T 1 

a(r') 

'AU gw (r') 

a(r) 

A log t' 

a(r') 

'AU g w (t') 


(r) L Alogr' J peak \k£ eak (T') 


= 2 + 


2 n - 16 


-1 r 


+ a - 


n — 2 
2n — 16 


AU gw (r) 

Alogr 


^peakW 


- peak 


n — 2 


-1 r 


AO gw (r) 
A logr 


peak 



(2n-16)/(n-2) 


i - 


(2n-16)/(n-2) 


-) 


2+(2n—16) / (n—2) 


1 - 


^peak( T form) ) 

/,A / \\ a-(2n-16)/(n-2)' 

/ K neakl T l \ 


(4.33) 


The scale factor is included in the integrant because the GW energy density decreases with 
time as oc a _1 (r) oc r -2 in the oscillation dominated era. We define r^ by the time the GW 
emission peak wavenumber reaches k, that is, &^ eak (rfc) = k. We implicitly assume r > r^. in 
this calculation, which means that Eq. (4.33) is valid for k > fc^ eak (r). 

The result of Eq. (4.33) at r = 50 is plotted as the black curves in Fig. 2 and well 
describes the GW spectra. Note that the GW energy density is roughly proportional to 

« (^feakA) (2n_16)/(n_2) for ^eak( T form) (= ^ e ak( T ) r / T form) > k > A^ eak (r) and n > 6. 
Thus, in principle, we can obtain the value of n from the observation of GWs with high 
frequencies. Note that for the case of n = 6 and 7, the power-law index of GW spectra is 
positive (16 — 2n > 0). This implies that the GW spectrum has a peak at r = Tf orm . This 
is because the emitted GW energy density of Eq. (4.31) decreases with time faster than the 
redshift oc a -1 for n = 6 and 7. This implies that the peak energy density emitted at r < Tf orm 
is smaller than the redshifted one emitted at r = Tf orm when the spatially averaged magnitude 
of the flat direction reaches the minimum of the potential. Thus, the GW peak wavenumber 
is determined by & P eak — ^p eak ( T form) ~ T form- Then the GW peak energy density decreases 
with oc a -1 due to the redshift until reheating completes. The redshift factor is given by 
®form/Q rh = (77RH/77form) 2 ^ 3 - This is the reason that the GW peak wavenumber is constant 
in time for the case of n = 6 in Fig. 2. The GW spectrum for n = 8 is too flat to determine its 
peak wavenumber precisely, though the GW spectra are well fitted by the above estimation. 

Next, let us consider the GW spectrum for very small wavenumber modes k < fc pea k- 
For such super-horizon modes, we expect that the leading contribution for the Fourier trans¬ 
formed transverse-traceless part of the anisotropic stress is independent of k due to the 
loss of causality at the large scale [14, 36]. As a result, the GW spectrum is proportional to 
k 2 for the modes smaller than but around fc peak , which is in agreement with our numerical 
simulations in Fig. 3. It is proportional to k for the modes much smaller than /c peak , which 
are out of reach of our numerical simulations. In addition, we have to take into account the 
effect of the reheating, at which the expansion law of the Universe changes. These effects on 
super-horizon modes leave their imprints on GW signals at present, which we discuss in the 
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next subsection. We can calculate the GW spectrum for super-horizon modes by substituting 
a constant TT t into Eqs. (B.6) and (B.7) in Appendix B and plot the results at f = 50 as 
the gray dotted curves in Fig. 2. 

From Fig. 2, we can see that GW spectra drop off for modes k > (S> fc pea k)> which 

enter the horizon before the spatially averaged magnitude of the flat direction reaches VEV 
(|</>|). This is because the produced energy density of GWs is proportional to the fourth power 
of the field value of the flat direction as shown in Eq. (4.28). The result depends nontrivially 
on ch osc because the dynamics of the flat direction around the origin of the potential is mainly 
determined by the Hubble induced mass term. Let us note that we are less interested in such 
high frequency modes in light of the detection of GW signals (see Fig. 5). 

The GW emission terminates around T de cay ; which is defined by |cjy osc | H 2 (t ) = m| (see 
Eq. (4.16)). This is because the flat direction starts to oscillate around the origin of the 
potential and cosmic strings disappear at that time. This implies that the GW spectrum 
has a peak (or at least bends) at the wavenumber corresponding to the soft mass of the flat 
direction. In fact, the GW peak energy density and frequency at that time is given by 

[^gw(r d ecay)] peak -14(2 + — —'j (\ 2(n _ 1} ^ ’ (4 ’ 34) 


^peak ^ g Wl(}) 
a('Tdecay) y/\ Ch osc 


(4.35) 


for the case of n > 8. Here, we determine the numerical prefactors from the results of our 
numerical simulations. Although GW spectrum has a peak at k ~ for the case of n = 6 
and 7, it bends at the wavenumber of Eq. (4.35) and has the energy density of Eq. (4.34). In 
either case, the GW spectrum contains information of the soft mass of the flat direction at 
that wavenumber because of the relation of Eq. (4.35). Hereafter, we denote the wavenumber 
defined by Eq. (4.35) as the peak wavenumber fc pea k- The shape of GW spectrum is in fact 
freezed out around r = T decay ~ 40, which is given by Eq. (4.16) with m^/Hi = 5 x 10 -4 and 
c HoBC = —15. After that, the peak wavenumber does not change and the GW energy density 
decreases as oc a -1 until reheating completes. 

Figure 4 is a schematic diagram for the resulting GW spectrum emitted from the cosmic 
strings. We calculate super-horizon modes, i.e., small wavenumber modes by using Eqs. (B.6) 
and (B.7) with a constant T(J t . The result is plotted as the blue dashed curve. The GW 
spectrum for small wavenumber modes near the peak wavenumber is proportional to k 2 , 
while that for the smaller wavenumber modes is proportional to k. We also take into account 
the effect of reheating, so that the GW spectrum bends at the wavenumber corresponding to 
the inverse of the horizon size at the end of reheating [17], 


,, Cn s - H(t Rn ). (4.36) 

a(tRH) 

We explain this phenomena in the next subsection in detail. We calculate large wavenumber 
modes by using Eq. (4.33) for the cases of n = 6, 8,10, and 12, and the results are plotted as 
the black curves. One can see that the GW spectrum for high frequency modes depends on 
the value of n. For n > 9, the GW spectrum has a peak at the wavenumber corresponding 
to the inverse of the horizon scale when cosmic strings disappear, that is, at the time of 
H(t) — m,p/^/\ ch osc \- For n = 6, 7, and 8, the GW spectrum does not have a peak, but bends 
at that wavenumber. Note that the magnitude of GW energy density depends on the value 
of the coupling in the superpotential A, though we use it to rescale the vertical axis in Fig. 4. 
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Log \k /{a H )]rescaled 

Figure 4. Schematic diagram for the GW spectrum emitted from cosmic strings. The horizontal 
and vertical axes are rescaled such that the GW energy density and the GW peak wavenumber are 
one. The black curves are plotted by using Eq. (4.33) for the cases of n = 6,8,10, and 12. The blue 
dashed curve is calculated by using Eqs. (B.6) and (B.7) with a constant TT t . We take into account 
a nonzero value of reheating temperature, which modifies the GW spectrum for wavenumber smaller 
than the inverse of the horizon scale at the end of reheating. 


4.4 GW signals at present 

The above calculations and discussions focus on the GW spectrum during the oscillation 
dominated era. In this subsection, we derive the GW spectrum at the present epoch, taking 
into account the redshift. 

The energy density of GWs at Td ecay , fl gw ( r decay) is given by our numerical simulations 
as shown and discussed in the previous two sections. Note that we consider the case that 
cosmic strings disappear before reheating completes, as we explained in Sec. 4.1. Since the 
GW energy density ft gw decreases with time as a -1 during the oscillation dominated era (i.e., 
the matter dominated era), its present value is given by 


Q ; ' 2 


“gw 


h (t 0 ) 


su 2 ( 

W(/rh)/ 

tu 2 ( JO^L) 

\g s (tRn)J 


4/3 


4/3 


g* (to) ) 
g*{tR hA 
g* {to) ) 


^gw(%H) 


Hru \ 

-^decay / 


2/3 


fig w (hdecay ): 


(4.37) 


where to is the present time, f l r h 2 (~ 4.15 x 10 -5 ) is the present energy density of radiation, 
and g* ( g s ) is the effective relativistic degrees of freedom for the energy (entropy) density. In 
the following calculations, we set g s {to) = 43/11, g*(to) = 2 + 21/11(4/11) 1 / 3 , and g s {tRH ) = 
5*(£rh) = 915/4 in the following calculations. Substituting [fl g w(7decay)] pca k i n t° Eq. (4.37), 
we obtain 


[figw/i 2 (to)] peak 


~ 3 x 10 -6 


( Icgpscl 1/2 
^ 10 3 GeV 


Trh 

10 9 GeV 


4/3 



2 n- 16 
n — 2 


-l 


( 1 

A 2 (n — 1) Mpj 


2/(n—2) 


(4.38) 
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where we use Eqs. (4.7), (4.8), and (4.34). The GW peak frequency at present /o is given by 


/o — 




9s {to) \ ^ / T) \ fcpeak 

9s{tmi)J \T RH ) 27ra(fRH) 

9s(to) \ 1/3 ( T 0 A / H rh A 2/3 fc peak 
5s(^Rh)/ \7rh/ \ Hdecay J 2?ra(t decay) 


Substituting Eq. (4.35) into this equation, we obtain 


fo ~ 700 Hz 


C H 0 , 


1—1/2 \ V 3 

I m 4> \ 


10 3 GeV 


Trh A 1/3 
10 9 GeV ) 


(4.39) 


(4.40) 


We should emphasize that this peak frequency depends only on |cfl- osc |~ 1//2 (~ H{T^ ecay )) 

and T rh and is independent of n and A. 

Here, let us consider the effect of reheating on the GW spectrum, which we mentioned 
in the previous section. In the absence of GW source, the metric perturbation is constant 
for super-horizon modes and decreases with a -1 for sub-horizon modes. Let us define aj~ 
as the scale factor at which GW with wavenumber k enters the horizon: k/ak = H{a,k)- 
For H{a) oc a ~ u , scales like a*, oc A; -1 /! 1 ' -1 ). Noting that H gw is proportional to the 

metric perturbation squared, we obtain H gw oc for sub-horizon modes, which 

have primordial GW spectrum of H gw before entering the horizon. If the primordial GW 
source is scale invariant, i.e., T? 7 = const, with respect to k, the primordial GW spectrum 
scales like Hgw oc k 5 . Since the expansion law of the Universe changes from H oc a~ 3 / 2 to 
H oc a -2 , we expect Hgw oc k and oc k 3 for GWs which enter the horizon before and after 
reheating, respectively, if the primordial GW source is scale invariant. As a result, the GW 
spectrum is expected to bend at the wavenumber corresponding to the horizon scale at the 
end of reheating [9, 17, 18]. Since we consider the case that cosmic strings disappear before 
reheating completes, the relevant GW with wavenumber corresponding to the horizon scale 
at the end of reheating is super-horizon mode when the GW emission terminates (r ~ Tdecay)- 
Thus, we can use T? t = const, in Eqs. (B.6) and (B.7) to estimate the GW spectrum due to 
the loss of causality for super-horizon scales [14, 36]. Then, using ji{x) 2 l l\x l / (2 1 + 1)! and 
ni(x) —> —{2l)\/{2 l l\x l+1 ) for x <C 1 in Eq. (B.8), we obtain H gw oc k for <C k <C T^ ay 
and H gw oc k 3 for k <C as expected. This means that the GW spectrum actually 
bends at the wavenumber around k — t^. We denote the frequency corresponding to that 
wavenumber at present as bend frequency /bend (see Eq. (4.36). It is calculated as 


, ( g.M \ 1/3 ( n \ 

/b “ d U((rh>1 vw 


^-bend 

2na{t R ii) 


~ 30 Hz 


Trh \ 
10 9 GeV ) 


(4.41) 


Note that this bend frequency is different from the one corresponding to k pe a k for n = 6 and 7. 
In these cases, GW spectrum bends at two frequencies: at / = /o (k = fc pea k) corresponding 
to the horizon scale at r = r^ecay and at / = /bend (k = ^bend) corresponding to the horizon 
scale at the end of reheating. 

Here, let us explain what kind of information we can obtain through the observation 
of GW spectrum emitted from cosmic strings. First, we can obtain the value of reheating 
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Figure 5. GW spectra generated by cosmic strings (dash and dot-dashed curves) and sensitivities 
of planned interferometric detectors (solid curves). We assume m$ = 10 2 GeV with Trh = 10 7 GeV 
(blue dashed curve), = 10 2 GeV with Xrh = 10 9 GeV (red dashed curve), and m$ = 10 3 GeV 
with Tr,h = 10 9 GeV (red dot-dashed curve). We also assume Ch obc = 10 and n = 10. We set the 
value of A so that the GW peak energy density is given by 0.01 at r = T^ecay- 


temperature through detection of the bend frequency of Eq. (4.41). Then we can obtain the 
soft mass scale of the flat direction through Eq. (4.40). We can also obtain the value of the 
power n because the slope of the GW spectrum for wavenumber larger than Eq. (4.40) is 
dependent on n (see Eq. (4.33) or Fig. 4). Finally, we can obtain the value of A through the 
measurement of the GW energy density Eq. (4.38). Thus, the observation of the whole GW 
spectrum gives us the values of |cjy osc | -1 / 2 m^, Trh, n, and |c # oac | -1//4 A . 

Figure 5 shows examples of GW spectrum predicted by the present mechanism. The 
GW peak energy density and peak frequency are given by Eqs. (4.38) and (4.40), which are 
obtained from our numerical results. While the shape of the GW spectrum for low frequency 
modes are calculated from Eqs. (B. 6 ) and (B.7) with a constant T-J t and tj = Tdecay; that 
for high frequncy modes are calculated from Eq. (4.33). Note that the high frecuency modes 
corresponding to the horizon scale at the time of cosmic string formation Tf orm are out of 
the range of Fig. 5, because the energy scale of inflation is much larger than the mass 
scale of the flat direction. We take = 10 2 GeV with Trh = 10' GeV (blue dashed 
curve), m 0 = 10 2 GeV with Trh = 10 9 GeV (red dashed curve), and = 10 3 GeV with 
Trh = 10 9 GeV (red dot-dashed curve). 

To discuss the detectability of GW signals, we plot single detector sensitivities for 
LISA [37] and Ultimate DECIGO [38] by using the online sensitivity curve generator in [39] 
with the parameters in Table 7 of Ref. [40]. Note that in some parameter space we can ob¬ 
tain Trh by using only one single detector for Ultimate DECIGO. One may wonder why our 
sensitivity curve is different from that in Ref. [18], where they claim that Ultimate DECIGO 
can reach well below D gw ~ 10 _lb for frequencies of 0.1 — 1Hz. This is because they assume 
taking two detectors correlation for a decade to obtain Trh by observation of the inflationary 
GW background, while we assume a single detector. A single detector is sufficient for our 
purpose because GW signals from cosmic strings are much stronger than inflationary ones. 
We also plot cross-correlation sensitivities for Advanced LIGO [41] and ET (ET-B configu- 
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ration) [42], assuming two detectors are co-aligned and coincident. In Fig. 5, we take the 
signal to noise ratio SNR = 5, the angular efficiency factor F = 2/5, the total observation 
time T = lyr, and the frequency resolution A/// = 0.1. CMB constraints (horizontal lines) 
are put on the integrated energy density of GWs f d log f£lgwh 2 (to) [43]. We find that GW 
signals would be observed by ET or Ultimate DECIGO. 

Note that we assume a very small A in Fig. 5 so that the future GW detector can observe 
GW signals. However, such a small A may result in a VEV (Eq. (4.4)) larger than the cutoff 
scale just after inflation, even though the VEV is smaller than the cutoff scale around the 
time of H(t) — m^,. In this case, the potential for the flat direction should be dominated 
by higher dimensional superpotential or Kahler potentials to stabilize its VEV at a smaller 
scale than the cutoff scale for a while. As the Hubble parameter decreases, only the lowest 
dimensional nonrenormalizable superpotential becomes relevant. Then the GW signals are 
determined by the discussion in this section. Finally, the potential is dominated by the soft 
mass term and the GW emission terminates. We should emphasize that the power-law index 
of GW spectrum for high frequency modes around the peak frequency depends on the value 
of n as shown in Fig. 4. Thus, for larger wavenumber modes, the GW spectrum bends at the 
wavenumber corresponding to the horizon scale when the potential becomes dominated by 
the lowest dimensional nonrenormalizable superpotential. However, such a high frequency 
modes are beyond the reach of future GW detecters, so that we neglect such effects in Fig. 5. 


5 Case B: W = 0 with {7(1) symmetry 


In this section, we consider the case that the potential for the flat direction is given by a higher 
dimensional term of Eq. (3.15), which comes from a nonrenormalizable Kahler potential. It 
is given by 


7nr(0) = a H H\t) 


101 


2m —2 


M 


2m —4 ’ 
PI 


where m > 3. The coefficient an is given by 


an — 


/MpA 

UJ 


2m—2 


(5.1) 


(5.2) 


where we assume a' H = 0(1). 

The following discussion and calculation in this section are essentially the same as the 
ones given in the previous section. In the next subsection, we briefly explain the dynamics 
of the flat direction in the early Universe. We numerically follow its dynamics and calculate 
its GW signal, and the results are shown in Sec. 5.2. Then we give physical interpretation of 
the results in Sec. 5.3. Finally, we discuss the detectability of GW signals in Sec. 5.4. 


5.1 Dynamics of the flat direction 

The whole potential for the flat direction V (</>) is given by 

I /i2m—2 

V(4) = ml H 2 + c H H\t ) H 2 + a H H 2 (t ) 1 -^^, (5.3) 

Mpi 

where we assume that ch is positive during inflation (c# inf > 0) and is negative ( ch osc < 0) 
after inflation. The positive c# inf makes the flat direction stay at the origin during inflation. 
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Then, after inflation ends, the negative ch obc makes it obtain a large VEV. The potential 
takes the minimum at 

/ | | \ l/(2m—4) 

M) = M P1 (~ M*), (5.4) 

\a H (m - 1 ) J 

as long as |c# osc | H 2 (t) m|. Here we have used |cff OBC | r^j (Mpi/M*) 2 (see Eq. (3.10)) and 

a// ~ (Mpi/M*) 2m ~ 2 (Eq. (3.16)). The VEV is as large as the cutoff scale M*. Note that the 
VEV of the flat direction is independent of the Hubble parameter and is constant in time. 
This is an important difference from the previous section (see Eq. (4.4)). 

Since the above potential possesses a global U( 1) symmetry, cosmic strings form after 
the end of inflation. By using numerical simulations (see the next subsection), we confirm 
that cosmic string network reaches a scaling regime well before cosmic strings disappear. 
In this regime, the number of cosmic strings in the Hubble volume is 0(1 ) and the energy 
density of cosmic strings p cs is given by Eq. (4.6). The energy density of cosmic strings per 
unit length p is given by Eq. (4.5). A typical width of cosmic strings is roughly given by 
the curvature of the potential, which is of the order of the Hubble length in the present 
case. This means that the width of cosmic strings increases with time. This is one of the 
outstanding characteristics of cosmic strings in this scenario compared with other cosmic 
strings considered in the literature. 

When the Hubble induced mass decreases down to the mass of the flat direction rn^. 
the flat direction starts to oscillate around the origin of the potential and cosmic strings 
disappear. The time when cosmic strings disappear is given by Eq. (4.7) as in the previous 
section. For the case around rxiA = 0(1) TeV, cosmic strings disappear much earlier than 
the BBN epoch, so that observational constraints on long-lived topological defects is not 
applicable to our scenario. As we explained in the previous section, we consider the well 
motivated case that the cosmic strings disappear before reheating completes. 

After they form at the end of inflation and before they disappear at t ~ fdecay: cosmic 
strings emit GWs. We numerically follow the evolution of the flat direction and calculated 
its GW signal. The results are shown in the next subsection. 

5.2 Results of numerical simulations 

In this subsection, we show our results of numerical simulations for the dynamics of the flat 
direction and its GW signal. 

By redefining the variables, we rewrite the equation of motion such as 

- w + (f"- - l) (?) * + (?) +(m ■ :> (?) |#,|2m_4 v =o<5 - 5) 

where 

(p = a^ 2vn Mp^a(t)(j) 

= (|c«J/(m- l)) 1/(2 ”- 4) a(t)0/(|«, (5.6) 

dr = Hidt/a(t ), (5-7) 

rh<!> = rri'p/Hi (< 1), (5.8) 

x = HiX. (5-9) 


- 24 - 




The parameter an can be absorbed by the redefinition of 4>, so that our numerical simulation 
can be applied to any case with an arbitrary value of an- 

We follow the time evolution of the flat direction by solving the equation of motion 
Eq. (5.5) numerically. Details of the numerical method are explained in the previous section. 
We confirm that each Hubble volume contains 0(1) cosmic strings and a typical width of 
cosmic strings is roughly proportional to the Hubble radius. From our numerical simulation, 
we find that the spatially averaged magnitude of cj) grows with time and eventually reaches 
of the order of the VEV (Eq. (5.4)) around tH, l = ?f orm — 10 for the case of ch osc = —15. 

We calculate the transverse-traceless part of the anisotropic stress from the numerical 
simulation and obtain the energy density of GWs using the equations in Appendix B. The 
results are plotted in Fig. 6. One can find that the GW peak energy density is almost constant 
in time and is almost independent of m. We give physical interpretation of these results in 
the next subsection. 

Taking a non-zero value of into account, we find that the flat direction starts to 
oscillate around the origin of the potential and cosmic strings disappear around r = r^ecay — 
40, which is given by Eq. (4.16) with m^/H* = 5 x 10 -4 and ch osc = —15. After that, the 
GW peak wavenumber becomes constant (blue dashed line in Fig. 6) and the GW spectrum 
decrease adiabatically as oc r^ 2 (oc a” 1 ). As we discussed in Sec. 4.2, the above choice of 
/ H, is convenient for numerical simulations, though it does not change the dynamics of flat 
direction until r ~ Td eC ay I n reality, it can be many orders of magnitude smaller depending 
on inflation models. 


5.3 Physical interpretations 

In this subsection, we discuss how we can understand the results of our numerical simulations 
shown in Fig. 6. 

Since during the scaling regime the typical length scale of cosmic string dynamics is the 
Hubble length, the GW emission peak wavenumber is typically given by the Hubble scale 
like Eq. (4.25). The produced energy density of GWs is estimated as Eq. (4.28) but with the 
VEV given by Eq. (5.4) in this case: 


' asv i /wy./ ic«j y /(m - 3> ( m . y 

A log rj peak V m pi/ \a H (m- 1)) \M P J 


(5.10) 


Cosmic strings emit GWs with a constant energy density expressed above. 

We plot our numerical results of AH gw /Alogr in Fig. 7. It shows that the produced 
energy density of GWs is constant in time as expected. We can determine the numerical 
prefactors for the emission peak wavenumber and the produced energy density of GWs from 
the numerical results like 


k A 

peak 

a(t) 


2.517(f), 


AH 


gw 


Alogr 


~ 10 


- peak 


\ c H n , 


a H (m - 1) 


2/(m—2) 


(5.11) 

(5.12) 


We find that the numerical prefactors are almost independent of m and r up to by a factor 
of two. 

Figure 7 shows that AH gw /Alogr is proportional to a certain power of k for large and 
small wavenumber modes like Eq. (4.32), where a — 2.5 in this case. The GW spectra for 
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Figure 6. GW spectra obtained from our numerical simulations. The energy density of GWs are 
rescaled by a H . The red, green, and magenta curves are GW spectra at the conformal time of 

T<nHi = 10,20, and 50, respectively. We also show the time evolution of the GW peak location as 
the blue solid (dashed) line for the case of m^/Hy = 0 (5 x 10 -4 ). We take m = 4 (upper left panel), 
6 (upper right panel), 8 (lower left panel), and 10 (lower right panel), fixing ch osc — —15. The black 
curves are plotted by Eq. (5.15) at TayHy = 50 and well describes large wavenumber modes. The gray 
dotted curves are estimated from the discussion of causality at TayHy = 50 and well describes small 
wavenumber modes. 


large wavenumber modes (k > kp eak (r)) in Fig. 6 can be understood by an integration like 
Eq. (4.33), where we can effectively take n —> oo in this case. Thus, we obtain 


^gw (r) 


1 



^ fcpeak(r)\ 2 

1 - 

( k y 

4 

A logr 

peak 

V k J 

\ ^peak(Tform) J 


+ (ot 


1 

'A figw(r)' 


( kpeak(^) \ 2 

1 - 

(k peak (T)\ a 2 


A logr 

peak 

V k J 

V k J 


where Ty. is defined by /c pe ak( T fc) = k. 


(5.13) 

(5.14) 

(5.15) 
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Figure 7. Derivatives for GW spectra in terms of the conformal time. The energy density of GWs 
are rescaled by the factor given by Eq. (5.10). We plot GW spectra for TUiHi = 20,30,40,50,60,80, 
and 100. We take m = 4 (upper left panel), 6 (upper right panel), 8 (lower left panel), and 10 (lower 
right panel), fixing ch obc = —15. 


The result of Eq. (5.15) at r = 50 is plotted as the black curves in Fig. 6. It is consistent 
with the large-wavenumber behaviour of GW spectrum found in our numerical simulations, 
though it is larger than our result by a factor of 2 in the case of m = 4. Note that the 
GW spectrum is expected to be proportional to k~ 2 for ^eak^form) > k » ^ eak (r). This 
asymptotic behavior is different from the one obtained in the previous section because the 
VEV is constant in time in the present case. The gray dotted curves in Fig. 6 show that 
simulated super-horizon modes are consistent with the spectrum expected from the discussion 
of causality [14, 36] as explained in the previous section. 

The blue dashed lines in Fig. 6 show the evolution of GW peak wavenumber in the case 
of ni'fi/Hi = 5 x 1CT 4 . Since the flat direction starts to oscillate around the origin due to its 
soft mass term, cosmic strings disappear at r ~ Td ecay . The GW peak energy density and 
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frequency at that time is given by 


If \ Ch I \ 2 /(™-2) 

[^gw(r d ecay)] peak - 10- 1} ) 

^peak ^ g ™(f) 

«(^decay) V\ C H osc \ 

where the numerical prefactors are determined from our numerical results. The shape of GW 
spectrum is in fact freezed out around r = r deca y — 40, which is given by Eq. (4.16) with 
m^/Hi = 5 x 10 -4 and ch osc = —15. After that, the peak wavenumber does not change and 
the GW energy density decreases as oc a -1 until reheating completes. 

In summary, the schematic form of resulting GW spectrum is similar to the one de¬ 
scribed in Fig. 4, except for the power-law index of large wavenumber modes. In the case 
of this section, the GW spectrum for large wavenumber modes is proportional to k~ 2 and is 
independent of m. This is because the potential minimum for the flat direction is constant 
in time. Thus, in principle, we can distinguish the cases with and without the higher dimen¬ 
sional superpotential through the observation of the GW spectrum with higher frequencies. 


(5.16) 

(5.17) 


5.4 GW signals at present 

The above calculations and discussions determine the GW spectrum during the oscillation 
dominated era. In this subsection, we derive the GW spectrum at the present epoch, taking 
into account the redshift. 

The peak energy density of GWs at r = T deca y, [n gw (T d ecay)] peak is given by Eq. (5.16) 
from our numerical simulations. Substituting [Hgw(^decay)] pcak into Eq. (4.37), we obtain the 
present value of GW peak energy density such as 


[Q gw /i 2 (to)] peak — 5 x 10 


( kgpJ 1/2 

l 10 3 GeV 


- 2/3 


Trh 


10 9 GeV 


4/3 


C H 0] 


a H (m - 1) 


2/(m—2) 

(5.18) 


where we use Eqs. (4.7), (4.8), and (5.16). From Eq. (4.39), the GW peak frequency at the 
present epoch is given by 


/o — 700 Hz 


c #oJ 1/2 m <j> \ / f t RH \ 1/3 
10 3 GeV j \10P GeV ) ' 


(5.19) 


Note that the peak frequency is independent of m and an- 

Since the GW spectrum is sensitive to what dominates the energy density of the Uni¬ 
verse, it contains information of reheating temperature as explained in the previous section. 
In fact, the GW spectrum bends at the frequency given by Eq. (4.41) corresponding to the 
wavenumber around k — r^ H [9, 17, 18]. 

Since there are three observables: Eqs. (4.41), (5.18) (see also Eq. (5.4)), and (5.19), in 
principle, we can obtain three parameters through observation of GW spectrum: \ch osc | _1//2 m^, 
Trh, and M*/Mpi (up to by a factor of order one). Again, let us stress that we can distin¬ 
guish the case of W ^ 0 in the previous section from the case of W = 0 in this section. This 
is because for high frequencies the GW spectrum scales as H gw oc yt 1 6-2n)/(n-2) ^ y- 2 ) for 
the former case and H gw oc /~ 2 for the latter case. 
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Figure 8. GW spectra generated by cosmic strings (dash and dot-dashed curves) and sensitivities 
of planned interferometric detectors (solid curves). We assume m$ = 10 2 GeV with Trh = 10 7 GeV 
(blue dashed curve), = 10 2 GeV with Xr H = 10 9 GeV (red dashed curve), and mn^ = 10 3 GeV 
with Trh = 10 9 GeV (red dot-dashed curve). We also assume ch obc = —15, a' H = 1, and m — 4. We 
set the value of M* as M 2 = M|j/10. 


If the cutoff scale was the Planck scale (M* ~ Mp\), we could obtain the mass of the flat 
direction and the reheating temperature Trh from observations of GWs (see Eqs. (4.41) 
and (5.19)). In this case, however, we can not ignore the effect of the backreaction of GW 
emission when we calculate the spectrum of GWs. In addition, the gradient energy density 
of the flat direction is comparable to the energy density of the Universe and thus it may spoil 
the isotropic evolution of the Universe. We can neglect these effects when the cutoff scale is 
less than the Planck scale. 

Figure 8 shows examples of GW spectrum predicted by the present mechanism. The 
GW peak energy density and peak frequency are given by Eqs. (5.18) and (5.19), which are 
obtained from our numerical results. While the shape of the GW spectrum for low frequency 
modes are calculated from Eqs. (B.6) and (B.7) with a constant T7 t and tj = Td eC ay> that En¬ 
large wavenumber modes are calculated from Eq. (5.15). In the figure, we take m$ = 10 2 GeV 
with T rh = 10 ' GeV (blue dashed curve), = 10 2 GeV with Trh = 10 9 GeV (red dashed 
curve), and rn^ = 10 3 GeV with Trh = 10 9 GeV (red dot-dashed curve). We could obtain 
rough scales of them from the observation of GW spectrum. To discuss the detectability of 
GW signals, we plot sensitivities of planned interferometric detectors. The sensitivity curves 
are the same as in Fig. 5, and details are given in Sec. 4.4. We find that GW signals would 
be observed by ET or Ultimate DECIGO. 


6 Case C: W = 0 without U( 1) symmetry 


In this section, we consider the same potential with the one considered in the previous section, 
Eq. (5.1), except for an additional 17(1) breaking term coming from higher dimensional Kahler 
potentials like 


C[/(i)b r 


- b H H\t) 
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( 6 . 1 ) 
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where m' > 2. We assume m! < 2m — 2 so that the potential for the flat direction is bounded 
from below. The coefficient bn is given by 


b H 


b\ 


H 



( 6 . 2 ) 


where we assume b' H = 0(1). Here we have redefined the phase of the field (j) such that b' H 
is postitive. This additional term explicitly breaks U( 1) symmetry, but there remains Z m i 
symmetry in that potential. 

The following discussion and calculation in this section are essentially the same as 
the ones explained in the previous two sections, though it is more complicated due to the 
additional term breaking 17(1) symmetry. In the next subsection, we briefly explain the 
dynamics of the flat direction in the early Universe. We numerically follow the evolution of 
its dynamics and its GW signal, and the results are shown in Sec. 6.2. Then we give physical 
interpretation of the results in Sec. 6.3. Finally, we discuss the detectability of GW signals 
in Sec. 6.4. 


6.1 Dynamics of the flat direction 

Throughout this section, we consider the following potential for the flat direction: 




2 m —2 


V(<t>) = m; j, \<f>\ z + c H H 2 (t) \cf>\ z + ^^(1)^=1 " ^H\t) 


0" 


pi 


Mi 


m '—2 
PI 


+ c.c. 


(6.3) 


where we assume that ch (= c H ]uf ) > 0 during inflation and ch (= ch osc ) < 0 after inflation. 
The positive c# ; f makes the flat direction stay at the origin during inflation. After inflation 
ends, the flat direction obtains a large VEV due to the negative ch osc - Note that the U( 1) 
breaking term unstabilizes the potential for the flat direction around the cutoff scale. For a 
small bn, the potential minimum is determined by the Hubble induced mass term as in the 
previous section (see Eq. (5.4)), while for a small c# inf , it is determined by the 17(1) breaking 
term. Thus, the VEV can be estimated as 


(0) = Ce 2p7ri/m, M P1 , 

"/ \ c „ I \ V (2m—4) / , b \ l/(2m-m'-2)' 

C ~ Max ( lCH ° sJ J 

\a H (m-l)J \m-la H J 
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~ M P1 ’ 


(6.4) 

(6.5) 

( 6 . 6 ) 


where p = 0,1,2,..., m' — 1. In the last line, we use Eqs. (3.10), (3.16), and (3.18). One can 
see that the VEV is as large as the cutoff scale M*. 

Since the nonzero VEV of the flat direction breaks Z m j symmetry, m! types of domain 
walls form after the end of inflation. In general, a typical width of domain walls is given by 
the inverse of the curvature of the potential at the minimum. In the present case, it is given 
by the inverse of the mass of the phase direction mg, which is given by 


m 2 e ~ (m') 2 b H H 2 (t)C m '~ 2 



Ml i 

Ml 



(6.7) 
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The energy of domain walls per unit area a can be estimated by 



~ (rn')“ i y^ : ^(t)C (m ' +2)/2 M| 1 (~ M*MpiH(t)). 


( 6 . 8 ) 

(6.9) 


These mean that the width of domain walls increases and the energy per unit area decreases 
with time. These are outstanding characteristics of domain walls in this scenario compared 
with other domain walls considered in the literature. By using numerical simulations (see 
the next subsection), we confirm that domain wall system soon reaches the scaling regime, 
in which the number of domain walls in the Hubble volume is of the order of the domain wall 
number m'. In this regime, the energy density of domain walls is roughly given by 

Pdw ~ m ' x a x H 2 x H 3 (~ . (6.10) 

We should emphasize that the energy density of domain walls never dominate that of the 
Universe. This is because the curvature of the potential is given by H(t) and decreases with 
time as oc a -3//2 , which is completely different from the ordinary domain walls [44]. 

When the Hubble parameter decreases down to the mass of the flat direction m^, the 
flat direction starts to oscillate around the origin of the potential and domain walls disappear. 
Since domain walls disappear well before the BBN epoch, our scenario is free from stringent 
constraints on long-lived domain walls [44], As we explained in Sec. 4.1, we consider the well 
motivated case that the domain walls disappear before reheating completes. 

After they form at the end of inflation and before they disappear at t ~ tdecayi domain 
walls emit GWs. In the subsequent subsections, we show our numerical results of domain 
wall dynamics and spectra of GWs emitted from these domain walls. 


6.2 Results of numerical simulations 

In this subsection, we show our results of numerical simulations for the dynamics of the flat 
direction and its GW signal. 

We redefine the variables in the same way as in the previous subsection (see Eqs. (5.6)- 
(5.9)) and obtain the equation of motion equivalent to Eq. (5.5) except for the additional 
1/(1) breaking term: 


ie"-VV+(w«-l)(S) ■?+(/) 

\ 4m—6 / /~ \ 2m'—2 

Vj M 2 m “V+ (^) m,_1 = o. 


( 6 . 11 ) 


The parameter an can be absorbed by the redefinition of <j) and bn, so that our numerical 
simulation depends on an and bn only though the combination of bn defined by 

bH^b H a^ m '^ 2)/{2m ~ 4) - (6.12) 


We follow the time evolution of the flat direction by solving the equation of motion 
Eq. (6.11) numerically. Details of the numerical method are written in Sec. 4.2. Figure 9 
shows one example of domain wall network evolution obtained by our numerical simulation. 
We set cn oac = —15, m = 4, and m! = 2. The green and cyan regions satisfy —7r/8 < 6 < 7r/8 
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Figure 9. Domain wall networks obtained by our numerical simulation at ra^IT = 20 (upper left 
panel), 50 (upper right panel), 80 (lower left panel), 100 (lower right panel). We set ch osc = —15, 
m = 4, bn = 5, and m! = 2. In the blue regions, the VEV of the flat direction satisfies \<j>\ < |(|</>|)| /30, 
where (|</f>|) is the VEV given by Eq. (6.6). In the green and cyan regions, the phase of the flat direction 
satisfies —7r/8 < 6 < n/8 and (7/8)n < 6 < (9/8)7r, respectively. The black line at the bottom-left 
corner represents the length of the comoving horizon size r at each time. 


and (7/8)-7t < 0 < (9/8)7r, respectively, where 6 is the phase of the flat direction. They 
represent the inside of domain walls. The blue region satisfies |0| < (|0|) /30 and describes 
cosmic strings in the limit of bn —> 0. The black line at the bottom-left corner in Fig. 9 
describes a unit of the comoving horizon size (given by r) at each time. One can see that a 
typical width of domain walls increases with time and is roughly proportional to the Hubble 
length. One can also see that each Hubble volume contains 0(1) domain walls of each type 
(m! = 2 in Fig. 9). From our numerical simulation, we find that the spatially averaged 
magnitude of cj) grows in time and eventually reaches of the order of the VEV (Eq. (6.6)) 
around Tf orm ~ 10 in this case. 

We calculate the transverse-traceless part of the anisotropic stress from the numerical 
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Figure 10. GW spectra obtained from our numerical simulations. The energy density of GWs 
are rescaled by a f/ 2,/ ' m 2 \ The red, green, and magenta curves are GW spectra at the time of 
TciiHi = 10, 20, and 50, respectively. We also show the time evolution of tire GW peak location as the 
blue solid (dashed) line for the case of rn^/Hi = 0 (5 x 10 -4 ). We take bn = 1 (upper left panel), 3 
(upper right panel), 5 (lower left panel), and 10 (lower right panel). The black curves at = 50 

are plotted by Eq. (5.15) and well describes large wavenumber modes. The gray dotted curves are 
estimated from the discussion of causality at raiHi = 50 and well describes small wavenumber modes. 


simulation and obtain the energy density of GWs using the equations in Appendix B. The 
results are plotted in Figs. 10 and 11. One can find that the GW peak energy density is almost 
constant in time. We give physical interpretation of these results in the next subsection. 

Taking a non-zero value of into account, we find that the flat direction starts to 
oscillate around the origin of the potential and domain walls disappear around r ~ id ecay , 
which is given by Eq. (4.7). After that, the GW peak wavenumber becomes constant and the 
GW spectrum decreases adiabatically as oc r -2 (oc a" 1 ). In Figs. 10 and 11, the GW spectrum 
is freezed out around r = Tdecay — 40, which is given by Eq. (4.16) with m^/Hi = 5 x 10^ 4 
and ch osc = —15. As we discussed in Sec. 4.2, the above choice of m^/Hi is convenient 
for numerical simulations, though it does not change the dynamics of flat direction until 
t ~ Tdecay- In reality, it can be many orders of magnitude smaller depending on inflation 
models. 
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Figure 11. Same as Fig. 10 but for the cases of m = 4 and m! = 3 (left panel), and m = 3 and 
to' = 2 (right panel). We assume bn = 5. 


6.3 Physical interpretations 

In this subsection, we discuss how we can understand the results of our numerical simulations 
shown in Figs. 10 and 11. 

During the scaling regime, the typical length scale of domain wall dynamics is the Hubble 
length. As a result, the emission peak wavenumber of GWs emitted from these domain walls 
is estimated by the Hubble scale like Eq. (4.25). Its quadrupole moment Q can be roughly 
given by [35] 


Q ~ H~ 2 x m'aH ~ 2 , (6.13) 

where m'oH ~ 2 is the total energy of domain walls within a Hubble volume. Such domain 
walls emit GWs with the luminosity of Eq. (4.27), which is proportional to the square of 
the quadrupole moment of domain walls. Thus, the produced energy density of GWs can be 
estimated as 


AH 


gw 


Alogr 


- peak 


b H c m ’ +2 



(6.14) 


Domain walls emit GWs with a constant energy density given by Eq. (6.14). Note that 
cosmic strings should coexist with domain walls because in the limit of bn —>• 0, the potential 
for the flat direction restores 17(1) symmetry. This implies that GWs are also emitted from 
cosmic strings. Thus, the total GW energy density is given by 


AH 


gw 


Alogr 


~ kh r b H C m +2 + . . 

peak \a H {m-l) 


C H 0 , 


2/(m—2) 


(6.15) 


where the second term comes from the contribution of cosmic strings and is relevant for 
small bn- We include a numerical prefactor k^ T to fit our numerical results. Note that the 
first term in the right-hand side of this equation is proportional to M^/Mpj (Eq. (6.14)), 
while the second term is proportional to Mf/Mp^ (Eq. (5.10)). This results from the fact 
that the energy density of domain walls is proportional to M* (see Eq. (6.10)), while that 
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of cosmic strings is proportional to M * (see Eqs. (4.5), (4.6) and (5.4)). Since the energy 
density of domain walls is much larger than that of cosmic strings, domain walls emit GWs 
more efficiently. 

We plot our numerical results of Afl gw /AlogT in Figs. 12 and 13. The results are 
constant in time as expected. We can determine the numerical prefactors for the emission 
peak wavenumber and the produced peak energy density of GWs from the numerical results 
like 


peak 

a(t) 


~ 2.5 H(t), 


An 


gw 


Alogr 


~ 10 

- peak 

fcbr — 0.4. 


k hr b H C m ' +2 + 


\ c H a 


a H {m - 1) 


2/(m—2) 


(6.16) 

(6.17) 

(6.18) 


We find that the numerical prefactors are independent of m, m/, and r, so that Eqs. (6.16) 
and (6.17) can be used for any values of m, m', and r. 

Figures 12 and 13 show that Afl gw /AlogT is proportional to a certain power of k for 
large and small wavenumber modes like Eq. (4.32), where a — 2.5 in this case. The GW 
spectra for large wavenumber modes in Figs. 10 and 11 can be understood by the same 
integration as Eq. (5.15). The result of Eq. (5.15) at r = 50 is plotted as the black curves in 
Figs. 10 and 11 and well describes the large wavenumber behaviour of GW spectrum. Note 
that the GW spectrum is proportional to k~ 2 for k fc pea k- 

The gray dotted curves in Figs. 10 and 11 show that simulated super-horizon modes are 
consistent with the spectrum expected from the discussion of causality [14, 36] as explained in 
Sec. 4.3. The blue dashed lines in Figs. 10 and 11 show the evolution of GW peak wavenumber 
in the case of m^/Hi = 5 x 10” 4 . Since the flat direction starts to oscillate around the 
potential minimum due to its soft mass term, domain walls disappear at r ~ Td eca y The 
GW peak energy density and frequency at that time is given by 


[n gw h (Tdecay)] 


peak 
^peak 


10 

T 


0.4^C m ' +2 + 


\ C H n 


a H {m - 1) 


2/(m—2) 


~ 3- 


m<t> 


a(Tdecay) \/|ctf osc | 


(6.19) 

( 6 . 20 ) 


where the numerical prefactors are determined from our numerical results. The shape of GW 
spectrum is in fact freezed out around r = T<j eC ay — 40, which is given by Eq. (4.16) with 
m^/Hi = 5 x 10” 4 and ch osc = —15. After that, the peak wavenumber does not change and 
the GW energy density decreases as oc a -1 until reheating completes. 

In summary, the schematic form of resulting GW spectrum is the same as the one de¬ 
scribed in the previous section, even though domain walls form. It is impossible to distinguish 
GW signals emitted from cosmic strings and domain walls. However, we can obtain infor¬ 
mation of SUSY scale and reheating temperature through the observation of GW signals, as 
we explain in the next subsection. 


6.4 GW signals at present 

The above calculations and discussions determine the GW spectrum during the oscillation 
dominated era. In this subsection, we derive the GW spectrum at the present epoch, taking 
into account the redshift. 
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Figure 12. Derivatives for GW spectra in terms of the conformal time. The energy density of GWs 
are rescaled by the factor given by Eq. (6.15) with a numerical constant k b r = 0.4, which we set to fit 
our results. We plot GW spectra for TUiHi = 20,30,40,50,60,80, and 100. We take bn = 1 (upper 
left panel), 3 (upper right panel), 5 (lower left panel), and 10 (lower right panel). 


The peak energy density of GWs at Td eC ay is given by Eq. (6.19) from our numerical 
simulations. Substituting [fi gw /i 2 (Tdecay)] pcak into Eq. (4.37), we obtain the present value of 
GW energy density such as 


[fl gw /r(t 0 )] peak — 5 x 10 


-7 / I c H a , 


1 — 1/2 \ -2 / 3 


10 3 GeV 
0 Ab H C m ' +2 + 


c H o: 


Trh \ 4/3 

10 9 GeV ) 


2/(m—2) 


a H (m - 1) 


( 6 . 21 ) 

( 6 . 22 ) 


where we use Eqs. (4.7), (4.8), and (6.19). From Eq. (4.39), the GW peak frequency at the 
present epoch is given by 


/o — 700 Hz 


cjtosd 1/2 m <t> \ / ( ^RH \ 1/3 

10 3 GeV J \10 9 GeV / ' 


(6.23) 
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Figure 13. Same as Fig. 12 but for the cases of to = 4 and m' = 3 (left panel), and to = 3 and 
to' = 2 (right panel). We assume bn = 5. 


Note that the peak frequency is independent of an, bn, m, and m!. 

Since the GW spectrum is sensitive to what dominates the energy density of the Uni¬ 
verse, it has information of reheating temperature as explained in the previous section. In 
fact, the GW spectrum bends at the frequency given by Eq. (4.41) corresponding to the 
wavenumber around k — r^ H [9, 17, 18]. 

There are three observables: Eqs. (4.41), (6.22), and (6.23), so that we can in principle 
obtain three parameters through observation of GW spectrum: |cff osc | _1,/2 Trh, and 
M*/Mpi (up to by a factor of 0( 1)). Let us stress that \ch osq \~ 1 ^ 2 and Trh are not 
affected by higher dimensional operators and are independent of the sources of GW signals, 
such as cosmic strings and domain walls. Note that we assume M* <C Mpi to avoid the effect 
of the backreaction of GW emission. 

Figure 14 shows examples of GW spectrum predicted by the present mechanism. The 
GW peak energy density and peak frequency are given by Eqs. (6.22) and (6.23), which are 
obtained from our numerical results. While the shape of the GW spectrum for low frequency 
modes are calculated from Eqs. (B.6) and (B.7) with a constant TL t and ry = r ^ ecay , that for 
high frequency modes are calculated from Eq. (5.15). In the figure, we take my = 10 2 GeV 
with T rh = 10' GeV (blue dashed curve), = 10 2 GeV with Trh = 10 9 GeV (red dashed 
curve), and my = 10 3 GeV with Trh = 10 9 GeV (red dot-dashed curve). To discuss the 
detectability of GW signals, we plot sensitivities of planned interferometric detectors. The 
sensitivity curves are the same as in Fig. 5, and details are given in Sec. 4.4. We find that 
GW signals would be observed by ET or Ultimate DECIGO. 

7 Summary and Conclusions 

We have investigated the dynamics of a flat direction in the early Universe, focusing on 
the case that the coefficient of the Hubble induced mass term is positive during inflation 
and is negative after inflation (oscillation dominated era). While the flat direction stays 
at the origin during inflation, it obtains a large VEV after inflation due to the negative 
Hubble induced mass term. In many cases, flat directions have 17(1) or Z n symmetries and 
thus topological defects, such as cosmic strings and domain walls, form just after the end of 
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Figure 14. GW spectra generated by domain walls (dash and dot-dashed curves) and sensitivities 
of planned interferometric detectors (solid curves). We assume m$ = 10 2 GeV with Trh = 10 7 GeV 
(blue dashed curve), = 10 2 GeV with T RH = 10 9 GeV (red dashed curve), and = 10 3 GeV 
with Trh = 10 9 GeV (red dot-dashed curve). We also assume c' H = 3/2, a' H = 1, m = 4, and m' = 2. 
We set the value of M* as M 2 = M-pJ 10. 


inflation. When the Hubble parameter decreases down to the soft mass of the flat direction, it 
starts to oscillate around the origin of its potential and topological defects disappear. These 
topological defects disappear well before the BBN epoch and they are free from constraints 
on long-lived topological defects. 

The topological defects emit GWs with an emission peak wavenumber corresponding 
to the Hubble scale at each time until they disappear when the Hubble parameter decreases 
down to the soft mass of the flat direction. This means that the resulting GW spectrum has 
a peak wavenumber corresponding to the soft mass of the flat direction [9]. Therefore, we 
can obtain the information of SUSY scale through the detection of GW signals. In addition, 
since the GW spectrum bends at the wavenumber corresponding to the horizon scale at the 
end of reheating, we can obtain the value of reheating temperature through the observation 
of the bending frequency. Since the GW energy density depends on the energy density 
of topological defects, it is related to the VEV of the flat direction. This means that the 
peak energy density of GW spectrum and the shape of GW spectrum for high frequency 
modes have information of higher dimensional terms for the flat direction. In this paper, 
we have calculated the GW spectrum emitted from cosmic strings and domain walls and 
have shown that the resulting GW spectra actually have the information of SUSY scale, 
reheating temperature, and higher dimensional terms for the flat direction. We investigate 
three cases for higher dimensional potentials for the flat direction; the one coming from a 
nonrenormalizable superpotential, the one coming from a t/(l) symmetric Kahler potential, 
and the one coming from 17(1) breaking Kahler potentials. Especially in the first case, we 
have found that the resulting GW spectrum for high frequnicy modes has information of 
power-law index for the nonrenormalizable superpotential. We have also found that the 
resulting GW signals would be measured by future GW observation experiments such as 
Advanced LIGO, ET, and Ultimate DECIGO. 

Here let us comment on the relation between the mass of the flat direction and Q- 
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ball formation [45-49]. Note that is a function of {(f)} through the renormalization. If 
dm^cj)} / dcj) < 0, the flat direction fragments into long-lived non-topological solitons called 
Q-balls soon after the flat direction starts to oscillate around <f> = 0 at t ~ tdecay- This is an 
obstacle in light of GW detection because it has been shown that GWs are diluted by the 
existence of Q-ball dominated era [50]. In this paper, therefore, we assume dm^cf) / dcf> > 0 
implicitly. This is realized in a wide class of gravity mediation models [48, 51]. In particular, 
Q-balls never form in SUSY models with a hierarchical spectrum between gauginos and the 
other sparticles in the MSSM. In models of gauge mediation, since the VEV of the flat 
direction induces the effective mass of gauge fields, soft masses are effectively suppressed as 
oc ( fr 2 for g \(f>\ > M s , where M s is a messenger scale [52, 53]. Thus, the flat direction need 
to contain a Higgs field to obtain a sizable value of and satisfy dm^cf)/dcf > 0 [54]. If 
GWs are observed in models of gauge mediation, we can identify as the Higgsino mass 
parameter fj,. 

Furthermore, let us comment on the baryon asymmetry (see footnote 5). The flat 
direction may be kicked by its A term along its phase direction after the inflation. The phase 
of the flat direction is randomly distributed through the phase transition. Therefore, though 
the baryon density may be locally generated by the Affleck-Dine mechanism, net baryon 
asymmetry is absent. The baryon asymmetry is thus never generated through the dynamics of 
this flat direction. To generate the baryon asymmetry, we may require another flat direction 
which has a large VEV during inflation and realizes the Affleck-Dine baryogenesis. The 
leptogenesis is also viable if the gravitino problem can be avoided [32-34, 55]. 

In the literature, it is known that cosmic strings emit Nambu-Goldstone modes effi¬ 
ciently, which leads to an effective friction on the dynamics of cosmic strings. As a result, 
cosmic string loops exist for a long time, which become more efficient source of GW emission 
than the cosmic string network in the ordinary cosmic strings [56-59]. On the other hand, 
in the present case, a typical width of cosmic strings is of the order of the Hubble length. 
Since the length of cosmic string loops can not be lower than its width, cosmic string loop 
formation is suppressed. This is why the effect of cosmic string loops can be neglected in the 
present paper. Note that the effect of cosmic string loops is incorporated in our numerical 
simulations. 

One might wonder if (tachyonic) preheating occurs for the flat direction just after the 
end of inflation and it affects the resulting GW signals [36, 60-65]. In fact, our numerical 
calculations incoperate the effect of tachyonic preheating for the fiat direction. However, let 
us emphasize that GWs in the observable band are not emitted from preheating which occurs 
when cosmic strings appear, but emitted from cosmic strings when they disappear. This is 
because the Hubble scale at the time of the preheating is much larger than the one at the 
time when topological defects disappear. Since GW signals with such high frequencies are 
beyond the detectability of future GW observations, we neglect GWs emitted from the pre¬ 
heating. We can also neglect the GW emission from the held / because its typical frequency 
is far beyond the observable region. This is because the characteristic energy scale of that 
dynamics is the mass of /, which is of the order of inflation scale. The typical GW frequency 
emitted from that dynamics is usually much higher than the observable region of 0.1-1 Hz. 
Furthermore, we consider that the reheating temperature is of the order of or less than 10 7 ~ 9 
GeV. This implies that inhaton decays into radiation perturbatively via Planck suppressed 
interactions. In this case, GW emission is not efficient at the time of reheating. 
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A Inflation models 

In this appendix, we briefly review two well-known inflation models in SUSY theories to 
clarify that the field whose F term drives inflation is generally different from the field whose 
oscillation energy dominates the energy density of the Universe after inflation. We use the 
Planck unit, that is, we set Mpi = 1 in this appendix. 

A.l Chaotic inflation model 

A simple chaotic inflation model in supergravity has been proposed in Ref. [6], where they 
introduced a shift symmetry to overcome the so-called r] problem in supergravity. They 
introduced two superfields X and I. The superfield I has a shift symmetry, which forbids 
the Hubble induced mass term for the imaginary part of its scalar component. This ensures 
the flatness of the potential for the imaginary part, so that it can be identified as inflaton. 
In order to introduce a potential for the inflaton, a small shift-symmetry breaking term is 
introduced in superpotential like 


W = mXI, 


(A.l) 


where X is another superfield and m is a parameter. 

In this appendix, we assume the simplest Kahler potential which satisfies the shift 
symmetry for the field I: 

K = - (I + I*) 2 + XX*. (A.2) 

This results in the following potential for the scalar components of the fields X and I: 


V(I, X) = m 2 e K 


\i \ 2 (i + |x| 4 ) +1 ac | 2 (i - \i \ 2 + (i + r ) 2 ( 1 + |I| 2 )) 


(A-3) 


Here, we can see that all scalar field except for the imaginary part of I are stabilized around 
the origin due to the exponential factor. On the other hand, the imaginary part of I, that is, 
the inflaton can take a field value larger than the Planck scale. Thus, the inflaton potential 
is simply given by 


V ~ JmV, (A.4) 

where x = Im//\/2 is the inflaton. In this way, a chaotic inflation model can be constructed 
in supergravity. Here, we should emphasize that the potential for inflaton Eq. (A. 4) does 
not come from the F term of the field / but from that of the field X. This results in a very 
important consequence explained in Sec. 3.2. This is also the case in variants of the above 
chaotic inflation model, e.g., the ones proposed in Refs. [66, 67]. 


- 40 - 




A.2 Hybrid inflation model 

In the simplest hybrid inflation model proposed in Ref. [25, 26], the superpotential is given 

by 


W = ASPP - n 2 S, 


(A.5) 


where S, P, and P are superfields, and A and /i are positive parameters satisfying A /i. 
The potential for their scalar components V(S, P, P) is thus given by 


V(S,*,V) 


JS |2 + |^| 2 + |^|2 


(1-|Sr + m|APP- M : 

2.T.* 


2|2 


+ \S\U A(1 + |P| Z )P-/x 2 P* + A(l + |P| 2 )P-/x 2 P 


.2,T <* 


(A-6) 
(A.7) 


When all scalar fields are well below the Planck scale (P, P, S <C 1) but satisfy the relation of 
S 3> 2^,/y/X, P and P are stabilized at the origin of the potential due to large effective masses. 
In this case, inflation is driven by a false vacuum energy density of ,u 4 . Driven by 1-loop 
induced potential for the field S , it rolls down to the origin of the potential. After the field 
S reaches the critical value of S cr = 2 / u/\/A, the water-fall fields P, P starts to roll down to 
the potential minimum at <[PP) = ^r 2 /A. Although this simplest model of hybrid inflation in 
supergravity might predict a spectral index inconsistent with the observation, many successful 
hybrid inflation models are considered in the literature based on this simplest model. Note 
that inflation is driven by the F term of the superfield S, while the energy density of the 
Universe is dominated by that of oscillation energy of water-fall fields P and \P after the end 
of inflation. In the main part of this paper, the field X is identified with the field S while I 
collectively represents water-fall fields \P and \P. 


B Method for the calculation of GWs. 

GWs are emitted by cosmic strings [10-13, 15, 71, 72] and/or domain walls [14-16] between 
the two phase transitions, that is, after the end of inflation and before the time of fdecay- la 
this appendix, we explain how to calculate the spectrum of GWs. The calculations are based 
on Refs. [14, 36], which is suitable to our situation compared with the method to calculate 
GW energy density from localized sources derived in Ref. [68]. 

The Fourier transformed transverse-traceless part of the spatial metric perturbation h t j 
obeys the following linearized Einstein equation: 

h'lj + + k 2 hij = 16v tGT? t , (B.l) 

where the prime denotes a derivative with respect to conformal time r defined by dr = a _1 dt, 
and is the Fourier transformed transverse-traceless part of the anisotropic stress. When 
the source term T? T lasts during the interval of [r;, Tf], the solution of Eq. (B.l) for Tf < t < 
t rh is given by [14] 

hij (r, k) = Aij (k) ^-ji (hr) + B {j (k) (fer), (B.2) 
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where j± and n\ are the spherical Bessel and Neumann functions of order one. The coefficients 
Aij and Bij are given by 


Aij (k) 


Bn (k) 


— I67 tG J dr' r / a(r / )ni(A:r / )T,^ T (rh k) , 
167 rG f dr r / a(r , )ji(fcr / )T^ T (r^k) . 


(B.3) 

(B.4) 


After reheating completes (r > trh), the solution is described by the spherical Bessel func¬ 
tions of order zero: 


hij (r, k) = A' t ■ (k) i— j 0 ( kr ) + B\ ■ (k) —n 0 {hr ). (B.5) 

J a a 

We match the solution of Eq. (B.2) with that of Eq. (B.5) at r = trh- We then obtain A’ 
and B l such as 


with 


Kj ( k ) = - 16 ttG J dr t’o,{t') (r'.k) , 
B^ (k) = 16v tG f d t T / a(r')/s(A:T / )Z;y T (r',k) , 


(B.6) 

(B.7) 


f A (kr') = [awi(kT') - a 2 ji(kr')] , 

f B {kT’) = [-bini(kT') + b 2 ji(kT')\ , (B.8) 

where the coefficients are given by 

a\ = x 2 [j 1 (x)d x n 0 (x) - n 0 (x)d x ji(x)\ \ x ^ k r Rii , 

a 2 = x 2 [ni{x)d x no(;x) - n 0 {;x)dxni(x)} |a^fcr RH , 

bi = - x 2 \ji(x)d x j 0 (x) - jo(x)d x ji(x)] \ x ^kr Rii , 

b 2 = - x 2 [n 1 (x)d x j 0 (x) - jo{x)d x ni(x)\ \ x ^ kTRH . (B.9) 

The energy density of GWs can be calculated from [35] 

Pgw = 22 (hij (^> x ) hij (t, x)) , (B.10) 

where (• • •) represents an ensemble average for a stochastic background. We replace the 
ensemble average by an average over a volume V of the comoving box. This is a good ap¬ 
proximation for sub-horizon modes. Using the following convention for the Fourier expansion: 

r (Uu 

M x ) = J J2^y$hij (k)e _ * k ' x , (B.ll) 

we rewrite the energy density of GWs as 


1 1 

Pgw “ 327rGa 4 V 


d 3 k k 2 
(2 tt) 3 T 


E(Ki 2 +i^r)- 


(B.12) 
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Here we have used h'- ~ khij and taken average over the oscillation period since we can 
observe only sub-horizon modes k -C r. Also we have averaged over time in Eq. (B.12) 
because we are not interested in the resolution of the oscillation of hij(r ) with time. We 
define the GW spectrum as 


H gw (r) 


1 dp gw (r) 
Ptot{r) dlog k ’ 


(B.13) 


where ptot(T~) (= 3Afpji/ 2 (r)) is the total energy density of the Universe. We can calculate 
the GW spectrum such as [36] 


k 5 

!igw(r) ~ 192rWIP 


d^E(KI 2 + |sij| 2 )- 


(B.14) 


We can calculate T? t from lattice simulation for the classical dynamics of a scalar field and 
then calculate the value of A' and B' through Eqs. (B.6) and (B.7). Finally, we can obtain 
the GW spectrum from Eq. (B.14). 


C Case with flavor symmetry 

In this appendix, we investigate the components of the flat direction in detail. In particular, 
we consider the case with flavor symmetry, where topological defects may not form when the 
flat direction obtains a large VEV and phase transition occurs. However, Nambu-Goldstone 
(NG) modes are randomly distributed over the Hubble scale through the phase transition, 
so that GWs are emitted from their relaxation dynamics within the Hubble horizon. The 
resulting GW spectrum is qualitatively the same as the ones obtained in the main part of this 
paper. This is a realization of so-called self-ordering scalar fields in SUSY theories [69-72]. 

In the next subsection, we investigate flavor components of the flat direction in detail 
and explain what happens in the case with flavor symmetry. Then we briefly review a method 
suitable to the calculation of GW spectra in this case, which has been derived in the literature 
in the context of self-ordering scalar fields. Using this method, we calculate the GW spectrum 
and show the results in Sec. C.3. 

C.l Flavor symmetry 

In order to investigate flavor components of flat directions, let us consider the LH U flat 
direction as an illustration. Since the left-handed lepton doublet L contains three flavors, 
the LH U flat direction has an ambiguity in terms of the flavor direction, that is, any flavor 
directions are D-flat directions. This implies that the LH U flat direction can be written as 


((L)\ \ 


/ «i 
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{L)\ 
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a 2 ■ 



{L)l 

“71 




V (Hu) 2 y 


\1 
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where the superscripts are SU(2) indices and the lower indices are flavor indices. Here, 
satisfy 

^h| 2 = L (C.2) 
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Any combinations of a* satisfying Eq. (C.2) are D-flat directions. This implies that the LH U 
flat direction has U( 3) symmetry if we take into account only D-flat condition. However, 
Hubble induced terms and higher dimensional superpotentials may break that symmetry as 
we can see in the following discussion. 

Next, let us consider Hubble induced mass for the LH U flat direction. If we assign R- 
charges as, for example, R{Q) = R(L) = R(e c ) = R(u c ) = R(d c ) = 3 and R(H U ) = R(H ( i) = 
—4, 11 nonrenormalizable Kahelr potentials take the following form: 


* = E iWii 2 + |ff»i 2 +£w 2 (i + E + % i« 


Jiu 

C X I TT 12 

M 2 l"«l 


+A ml „(/) (l+E W (L) ' (L)i + |H ' 


M 2 


(C.3) 


where Ca (a = A, /) are Hermitian constant matrices and c^ u (a = A, /) are real constants. 
As explained in Sec. 3.2 and Appendix A, the minimal Kahler potential for the field I is given 
by Amin (I) = (/ +1*) 2 / 2 for the case of the chaotic inflation model proposed in Ref. [6] and 
A m in(-0 = \I\ 2 for the case of hybrid inflation models. The F term of A drives inflation and 
satisfies the relation of \8W/dX\ 2 ~ 3H 2 ni Mp V The above Kahler potential induces Hubble 
induced mass terms for fields tpi such as 


= L 4 h ,FFm)i +eg- h 1 i h, 

i 

M 2 


2 

u\ ? 


^ — 
-Hi nf 


= 3 1 - c 




X 


M 2 


-'Anf 


, rr Mp, 

— Q (1 _ _ £1 

“ 1 x M 2 


(C.4) 

(C.5) 

(C.6) 


during inflation. After inflation ends, the energy density of the Universe is dominated by the 
oscillation of the field I, so that Hubble induced mass terms are given by 


vh ~ = E 4„,Au>, , (A + i//„i 2 , (c.7) 

i 



The coefficients c a ’s depend on underlying high energy models beyond the cutoff scale M*. 
Without any symmetry, we expect that c a ’s are not related with each other. This means that 
Hubble induced mass terms usually break the flavor symmetry, which is assumed in the main 
part of this paper (see footnote 4). However, one may assume flavor-universal coefficients 
such as c l a = c^<W , where Sij is Kronecker’s delta, just like the case in the constrained MSSM, 
for instance. In this case, the LH U flat direction has U( 3) symmetry. 

Nonrenomalizable terms should stabilize the flat direction as we discussed in Sec 3.3. It 
can take a form of 


Anr = a 


H 


l H 2 {L)* il {L) jl ... (L)* m _ n (L) jrn _ n |A, 


«i 2n ~ 2 /M, 


2m—4 
PI J 


11 See footnote 6 for discussion about effects of R-symmetry breaking. 
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where the coefficient an is of the order of (Mpi/Af*) 2m-2 and may also break flavor symmetry 
in general. If we consider, for example, the case of a^' -* m —— aH §*di .. tfm-njm-n^ 12 
the superpotential respects 17(3). 

Suppose that the Hubble induced mass term is flavor-universal and the higher dimen¬ 
sional superpotential has the 17(3) symmetry. As we can see from Eqs. (C.5) and (C.8), 
the coefficient of the Hubble induced mass term during inflation, c# inf , is generally different 
from the one during the oscillation dominated era, c# oac [9]. We consider the case that the 
coefficient of the Hubble induced mass term is positive during inflation and is negative after 
inflation. In this case, the flat direction stays at the origin during inflation and has a large 
VEV after inflation, that is, the phase transition occur just after inflation. At this phase 
transition, the 17(3) symmetry is broken to 17(2) symmetry. In this case, topological defects 
do not form at this phase transition. However, there is another source of GWs. Since NG 
models have no correlation over distances larger than the Hubble radius due to the causality, 
they are randomly distributed over the Hubble scale. Therefore, the NG modes obtain gradi¬ 
ent energy. The gradient energy turns into the kinetic energy and decreases inside the Hubble 
radius. As a result, the modes with wavenumber around the Hubble scale (k ~ aH(t )) domi¬ 
nantly emits GWs [69-72], Since the NG modes outside the Hubble radius grow through the 
nonlinear couplings, the peak frequency of the resultant GW signals correspond to the Hub¬ 
ble scale when the flat direction start to oscillate around the origin at f ~ idecay (Eq. (4.7)). 
This is another source of GWs which gives a contribution to the GW signal. Note that it is 
possible that, say, there is 0(N ) x U (1) symmetry and is broken down to 0(N — 1) symmetry. 
In this case, topological defects also form and the resulting GW signal is given by the sum 
of both contributions. In the following subsections, we calculate the GW energy density and 
spectrum from the above mechanism. 

Hereafter we focus only on the global symmetry of the flat direction for the following 
reasons. The gauge fields may be randomly distributed at the phase transition. We can take 
the unitary gauge if the topological defects are not formed. Since, the broken gauge field 
obtains the effective mass of g (q i), the field configuration begins to oscillate and becomes 
trivial ( A^ = 0) soon after the phase transition. For the unbroken gauge symmetry, the 
corresponding gauge field remains massless. Here we should note that while long wavelength 
modes of such gauge fields remain constant until they enter the horizon, the NG modes grow 
through the nonlinear term [72]. Therefore, GW emission from the local symmetry can be 
neglected. 

C.2 Analytical calculation of GW spectra for self-ordering scalar fields 

Here we explain a method suitable to the calculation of GW spectra emitted from self-ordering 
scalar fields. We consider the case with 0(N ) symmetry for simplicity. Note that LH U flat 
direction considered in the previous section has 0(6) symmetry rather than 17(3) symmetry 
if gauge and Yukawa couplings are ignored. This is similar to 0(4) custodial symmetry in 
the Standard Model. We expect that the results do not change qualitatively for the case 
with other symmetries because the discussion in the previous subsection hardly depends on 
the detail of the manifold where NG bosons reside. 

Instead of replacing the ensemble average by an average over a volume of the comoving 


12 This assumption is inconsistent with the observations of neutrino oscillations [73]. In this appendix, 
however, we take the LH U flat direction as an example to take advantage of its simple flavor structure. 
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box in Eq. (B.10), here we just rewrite it by assuming 


(^ T (n, k)T? T *(T 2 ,k')) = (2nfd 3 (k-k / )U(k,T 1 ,T 2 ), (C.ll) 

and then obtain the energy density spectrum of GWs in the radiation dominated era after 
reheating as 


2-jik 5 f Tt f T{ 

WT)= 3WWMi,j, dT1T1 “ (Tl) /, d75T2 “ W 

x [fA{kTi)f A (kT 2 ) + fB{kTi)f B (kT 2 )] n(fc,n,r 2 ) . (C.12) 


To calculate the spectrum of GWs, we quote the result of n(A;, 71,72) calculated in the 
literature [69, 71, 72], where they have considered a global 0(N ) symmetric scalar field with 
a quartic potential. While we consider higher power potentials for flat directions, their result 
is applicable to the present case. This is because the potential for scalar fields is irrelevant 
to NG bosons after the scalar fields obtains VEV. Within the large N approximation, they 
have obtained 


n(fc,r 1 ,r 2 ) 


V ab (k,T 1 ,T 2 ) 



l-(k-q) 2 


x 


V ab {q, 71, T 2 )V ab {\k - q| , 7 - 1 , 75 ), 
2835 3 , 3 / 2 J3(kTi)J 3 (kT 2 ) 5 ab 

16 12 ( kTi) 3 (kr 2 ) 3 N 



(C.13) 


(C.14) 


where Js(x) is the Bessel function of order 3. 

C.3 Results for self-ordering scalar fields 

When the Hubble parameter decreases down to the mass of the flat direction m^, the flat 
direction starts to oscillate around the origin of the potential. After that, the field value of 
the flat direction decreases with time due to the expansion of the Universe such as <f> oc a -3//2 . 
After the phase transition at the end of inflation and until they start to oscillate around the 
origin at t ~ fdecay, the dynamics of NG modes emit GWs continuously. We obtain the GW 
spectrum by evaluating Eqs. (C.12), (C.13), and (C.14) numerically. 

We obtain GW peak energy density and peak wavenumber at i ~ fdecay such as 

flgwh 2 (Tdecay) ^ 0.141V -1 > 

^peak ~ 2 3 
fl( 7 "decay ) \/|c,ff oac | 

Note that we assume M* /M-p\ < 1 so that we can neglect the backreaction of the flat direction 
dynamics on the cosmic expansion. Substituting the above formulas into Eqs. (4.37) and 
(4.39), we obtain the present value of the GW peak energy density such as 


(C.15) 

(C.16) 


tt gw h 2 (to) ~ 3 x 10 8 N 1 


( \ c HoJ 1/2 m# 
^ 10 3 GeV 


Trh 

10 9 GeV 
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and the present value of peak frequency /o such as 


Trh 

10 9 GeV 


1/3 


fo ~ 550 Hz 


/1 1-1/2 \ 3 / 3 

( \ c HoJ 

^ 10 3 GeV y 


(C.18) 


Since a typical length scale of the dynamics of NG modes is the Hubble length, the 
emission peak wavenumber of GWs is estimated as the Hubble scale. In addition, the GW 
energy density of Eq. (C.12) is proportional to (0) 4 . These observations imply that the 
resulting GW spectrum emitted from the dynamics of NG modes is qualitatively the same as 
the ones obtained in Secs. 4, 5, and 6, depending on higher dimensional potentials for the flat 
direction. In particular, since the GW spectrum is sensitive to what dominates the energy 
density of the Universe, it has information of reheating temperature as explained in those 
sections. In fact, the GW spectrum bends at the frequency given by Eq. (4.41) corresponding 
to the wavenumber around k — [9, 17, 18]. 

In summary, the schematic form of resulting GW spectrum is similar to the one described 
in the main part of this paper, even if there is a large flavor symmetry such as O(N). 
Figure 15 shows examples of GW spectrum generated by the self-ordering scalar fields. We 
assume that (< f>) = M* is constant until t = fdecay as in Secs. 5. The GW spectra are 
calculated from Eqs. (C.12), (C.13), and (C.14). In the figure, we take rn^ = 10 2 GeV with 
Trh = 10 7 GeV (blue dashed curve), = 10 2 GeV with Trh = 10 9 GeV (red dashed 
curve), and m^ = 10 3 GeV with Trh = 10 9 GeV (red dot-dashed curve). To discuss the 
detectability of GW signals, we plot sensitivities of planned interferometric detectors. The 
sensitivity curves are the same as in Fig. 5, and details are given in Sec. 4.4. Detecting these 
GW signals seems somewhat challenging in the current design of Advanced LIGO, ET, and 
Ultimate DECIGO unlike those from topological defects in Secs. 4, 5, and 6. However, let 
us stress that if we can detect these GW signals from somehow, for example, through longer 
observation time T > 1 yr, they provide us invaluable information of |c/7 osc |~ l/2 m^, Trh, 
and (up to by a factor of 0(1)). 
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